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ABSTRACT 

In  salvo  firing,  the  smallest  ballistic  dispersion  is 
not  always  most  desirable.   Deliberate  increase  of  the 
ballistic  dispersion  can  improve  the  probability  of  des- 
troying the  target.   Our  concern  in  this  dissertation  is 
the  optimization  of  such  "artificial11  dispersion  in  two- 
dimensional  salvo  models.   In  some  cases  no  closed  form 
solution  is  available,  but  we  are  able  to  offer  efficient 
methods  for  the  computation  or  approximation  of  the  salvo 
kill  probability.   In  other  cases  we  are  able  to  derive 
approximate  formulae  for  the  optimal  ballistic  dispersion. 
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I.   INTRODUCTION 

Let  us  begin  with  a  classical  example  presented  in  Chap- 
ter 6  of  the  textbook  by  Morse  and  Kimball  [Ref .  1] .   Suppose 
an  airplane  carries  two  bombs  to  attack  a  railroad  track. 
The  plane  flies  along  a  course  perpendicular  to  the  track, 
and  drops  the  bombs.   There  are  three  ways  to  drop  two 
bombs:   (1)  together  in  salvo,  (2)  spaced  a  suitable  distance 
apart  aiming  the  midpoint  of  the  pattern  at  the  center  of 
the  track,  and  (3)  dropping  each  on  a  separate  run  over  the 
target. 

We  intend  to  illustrate  the  method  of  calculating  the 
probability  of  track  destruction  in  the  three  cases,  the 
purpose  being  to  determine  which  method  of  attack  is  the 
best. 

Let  us  consider  Case  (1) ,  salvo  bombing.   The  two  bombs 
leaving  the  plane  simultaneously  will  hit  the  ground  at  some 
distance  apart;  the  impact  point  being  random.   For  the  time 
being  the  random  vector  of  the  impact  position  from  the 
center  of  impact  is  termed  the  ballistic  error.   The  aiming 
is  also  not  free  from  error.   The  plane  aims  at  a  point  on 
the  center  of  the  track,  but  the  center  of  impact  will 
deviate  from  the  aim  point.   This  deviation,  or  aiming  error, 
is  also  random.   The  impact  point  of  each  bomb  is,  there- 
fore, composed  of  two  errors;  first  the  aiming  error  which 


is  common  to  both  bombs,  and  second  the  ballistic  error 
which  varies  from  one  bomb  to  another. 

If  the  variance  of  the  ballistic  error  is  small  while 
that  of  the  aiming  error  is  large,  the  two  bombs  will  land 
almost  at  the  same  point.   Therefore  it  is  expected  that 
either  both  will  hit  or  neither  will  hit  the  track.   To  avoid 
the  latter  situation  and  to  improve  the  probability  of  des- 
truction, it  is  better  to  spread  the  landing  pattern  of  the 
bombs.   One  way  is  the  method  of  pattern  bombing  described 
as  Case  (2)  ,  another  is  salvo  bombing  employing  bombs 
with  a  suitable  value  of  variance  of  the  ballistic  error. 
Our  main  concern  in  this  dissertation  is  maximization  of  the 
probability  of  salvo  destruction  by  choosing  a  suitable 
value  of  the  variance  of  the  ballistic  error.   In  the  follow- 
ing we  call  it  the  optimal  ballistic  dispersion. 

The  problem  described  above  is  not  a  new  one.   In  an 
interesting  talk  before  a  conference  at  the  Ballistic  Re- 
search Laboratories  held  in  March  19  55  iRef .  2] ,  Merritt 
and  King  stated  that  an  approximate  formula  for  the  optimal 
ballistic  dispersion  was  derived  by  two  Englishmen  as  early 
as  in  19  36.   Following  the  work  in  WW  II,  a  number  of  arti- 
cles  were  published  on  both  the  salvo  and  pattern  firing 
models,  but  since  the  early  19  70's,  it  seems  that  these 
models  have  attracted  less  attention.   Further  detail  can  be 
found  in  two  excellent  review  papers  by  Eckler  [Ref .  3]    and 
Eckler  and  Burr  [Ref.  4] .   However,  little  is  known  as  to  the 


optimization  of  the  ballistic  dispersion  in  salvo  models. 
It  is  the  calculation  or  approximation  of  the  optimal 
ballistic  dispersion  that  is  the  subject  of  this  dissertation. 

A.   SALVO  FIRING 

In  the  following  chapters  we  deal  with  only  two-dimensional 
salvo  models.   The  reason  we  adopt  the  two-dimensional  model 
is  that  it  is  the  most  frequent  case  and  plays  an  important 
role  in  real  world  applications.   One-  or  three-dimensional 
theory  can  also  be  developed  to  parallel  our  investigation. 

Suppose  that  there  is  a  target  in  a  two-dimensional  space, 
and  a  salvo  of  n  weapons  is  delivered  against  the  target. 
Delivery  error  relative  to  the  target  is  assumed  to  be  com- 
posed of  two  parts,  the  aiming  error  and  the  ballistic  error. 

First  we  aim  at  the  target.   With  this  aiming,  the  center 
of  impact  point  of  n  weapons  is  determined.   Let  us  adopt  a 
Cartesian  coordinate  system  lx,yl  such  that  the  origin  coin- 
cides with  the  center  of  impact.   The  aiming  cannot  be  per- 
fect.  Let  the  components  of  the  aiming  error  be  -U  and  -V. 
Then  the  position  of  the  target  with  respect  to  our  coordinate 
system  is  given  by  (U,V) .   It  is  assumed  that  the  joint 
probability  density  of  U  and  V  exists  and  is  given  by  f..  (u,v)  . 

Now,  a  salvo  of  n  weapons  is  fired  against  the  target 
after  (U,V)  takes  a  value  (u,v)  which  is  unknown  to  us  and 
only  is  predictable  in  a  probabilistic  sense.   The  components 
of  the  impact  point  of  the  ith  weapon  are  denoted  as  X.  and 
Y..   It  is  assumed  that  (X.,Y.)  are  independent  and  identically 


distributed  random  variables  with  joint  density  f«(x,y) . 
In  Fig.  1,  a  typical  geometry  is  shown.   The  target  is  at 
(u,v) ,  and  the  impact  points  of  four  weapons  are  scattered 
around  the  origin. 

As  to  the  target  we  adopt  the  so-called  point  target 
concept:   It  is  assumed  that  the  target  is  either  completely 
destroyed  or  else  undamaged  by  each  weapon.   We  neglect  any 
possible  partial  damage  and  its  cumulative  effect.   The 
probability  that  a  weapon  landing  at  (x,y)  destroys  the  tar- 
get at  (u,v)  is  a  function  of  u-x  and  v-y,  denoted  as 
D(u-x,v-y)  and  is  called  the  damage  function. 

The  probability  of  destroying  the  target  is  calculated 
in  two  steps.   First,  suppose  that  the  random  variables 
U  and  V  take  some  values  u  and  v.   The  conditional  proba- 
bility that  the  ith  weapon  destroys  the  target  given  U  =  v 
and  V  =  v  is 


oo      oo 


p(u,v)   =   /    /    D(u-x, v-y) f2 (x,y) dxdy  .  (1.1) 


— OO     —CO 


The  impact  points  X.  and  Y.  are  all  assumed  to  be  independent, 
so  the  conditional  survival  probability  of  the  target  given 
U  =  u  and  V  =  v  is  {1  -  p(u,v) }  .   Therefore,  we  obtain  the 
probability  that  the  target  is  destroyed  by  the  salvo  of  n 
weapons — salvo  kill  probability — by  averaging  the  conditional 
survival  probability  with  respect  to  the  distribution  of  u 


0      (x4,y4) 


o     (x3^y3) 


0       (x1,y1) 


+       (u,v) 


0       (x2,y2) 


X 


Figure   1 
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and  v,  and  then  subtracting  it  from  unity: 


CO       CO 

n 


p=l-J    /   {1  -  p(u,v)}  f1(u/v)dudv  .  (1.2) 


—  CO      —CO 


B.   UPPER  BOUND  ON  THE  SALVO  KILL  PROBABILITY 

In  the  example  mentioned  earlier,  Morse  and  Kimball  ob- 
served that  the  probability  of  destroying  the  track  by  drop- 
ping two  bombs  in  salvo  is  always  smaller  than  the  destruc- 
tion probability  of  Case  (3) ,  namely  the  case  in  which 
the  two  bombs  are  dropped  independently  on  separate  passes. 
Similar  observations  are  also  pointed  out  in  the  review  paper 
by  Eckler  and  Burr.   The  observation  is  valuable,  but  there 
has  been  no  proof  given  to  this  fundamental  property  dis- 
covered numerically.   Indeed,  we  have  the  following  proposi- 
tion without  any  specific  assumptions  on  f, (u,v) ,  f~(x,y) 
and  D(x,y) . 

Proposition  1.1 

Let  the  salvo  kill  probability  with  n  weapons  be  denoted 
as  Pn 


CO       CO 


Pn   =   1  -  /    /   {1  -  p(u,v)  }nf1(u/v)dudv  .  (1.3) 


—CO     —CO 


Then, 
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Pn   <   1  -  (1  -Px)n  .  (1.4) 


Namely,  salvo  P   is  always  smaller  than  the  kill  probability 
of  n  independent  tries. 

Proof 

The  proof  is  immediate  from  a  theorem  in  the  book  by 
Hardy  and  others  [Ref .  5] .   Define  a  mean  M  of  a  non-negative 


function  g(x)  for  a  positive  r  as, 


Mr(g)   =  (  /  grf  dx} 


1/r 


where  f  =  f (x)  is  a  weighting  function,  positive  everywhere 
and, 

/  f  dx   =   1  . 

Theorem  19  2  from  Ref.  5  tells  us  that  if  0  <  r  <  s  and 
M  (g)  is  finite,  then 

9 


M    <   M   , 
r      s 


unless  g(x)  is  a  constant. 

The  weighting  function  f  corresponds  to  f, (u,vl  in  our 
proposition,  and  g  corresponds  to  our  l-p(u,v).   Therefore, 
for  n  >  1 ,  we  have 
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M1(l-p)  <    Mn(l-p) 


or 


(M^l-p)}11   <   {Mn(l-p)}: 


which  means 


(1-P,)n   <   1  -  P.  G 

1  n 


Here  is  an  outline  of  the  succeeding  chapters:   In  Chap- 
ter II,  we  investigate  the  simplest  salvo  model  which  is 
characterized  by  circular  normal  errors  and  circular  Gaussian 
damage  function.   After  that,  an  approximate  formula  for  the 
optimal  ballistic  dispersion  is  introduced  in  Chapter  III; 
we  regard  the  material  in  Chapter  III  as  the  "core"  of  this 
dissertation.   In  Chapter  IV,  the  approximate  formula  is 
applied  to  the  salvo  model  of  Chapter  II,  and  its  accuracy 
is  studied.   In  Chapter  V,  an  approximate  method  of  calcu- 
lating the  salvo  kill  probability  of  the  salvo  model  with 
the  general  damage  function  is  presented;  this  case  is  more 
difficult  than  the  one  where  the  damage  function  is  circular 
normal,  so  we  resort  to  further  approximations.   In  Chapter 
VI,  we  investigate  the  salvo  model  with  the  so-called  cookie 
cutter  damage  function,  relying  heavily  on  the  results  of 
Chapters  III  and  V.   Finally,  salvo  model  without  circular 
symmetry  is  discussed  in  Chapter  VII.   Our  concern  throughout 
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is  the  optimization  of  artificial  dispersion.   In  some  cases 
no  closed  form  solution  is  available,  so  we  offer  efficient 
methods  for  the  computation  or  approximation  of  salvo  kill 
probability  with  a  given  dispersion.   In  other  cases  we  are 
able  to  offer  closed  form  solutions. 
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II.   SIMPLE  SALVO  MODEL 

In  this  chapter  we  deal  with  the  simplest  salvo  model 
with  circular  normal  errors  and  a  circularly  symmetric 
Gaussian  damage  function. 

Let  us  adopt  a  Cartesian  coordinate  system  (x,y)  the 
origin  of  which  coincides  with  the  center  of  impact  of  the 
n  weapons.   The  target  position  with  respect  to  this  coor- 
dinate system  is  denoted  as  U  and  V  in  the  x  and  y  direction, 
respectively.   (U,V)  is  a  random  bias  common  to  all  the  n 
weapons.   In  this  chapter  (U,V)  is  assumed  as  a  circular  normal 
variate;  namely,  we  assume  that  U  and  V  are  independent  and 

identically  distributed  normal  random  variables  with  mean  0 

2 
and  variance  a,.   The  joint  density  of  U  and  V  is,  therefore 

2  ,  2 
u  +v 


The 


1         2°1 

f,  (u,v)   =   J  e        •  (2.1) 

27TQ, 

impact  point  of  the  ith  weapon  is  denoted  as  (X.,Y.) . 


It  is  assumed  that  (X.,Y.),  i  =  l,2,...,n  are  independent  and 
identically  distributed  circular  normal  random  variables 
with  joint  density  function 


2^  2 
x  +y 

2~ 

1       2°2 
f2(x,y)   =   -^  e      *       .  (2.2) 

2tto2 
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As  to  the  damage  function,  we  assume  that  the  so-called 

Gaussian  damage  function  with  circular  symmetry:  If  the 

target  is  at  (u,v)  and  a  weapon  impacts  at  (x,y) ,  then  the 
destruction  probability  is  given  by 


2       2 
(u-x)  +(v-y) 

:  2 

D(u-x,v-y)   =   e        za  .  (2.3) 


In  the  following,  let  us  call  this  model  the  Simple  Salvo 
Model. 

According  to  Grubbs  [Ref.  6],  the  Simple  Salvo  Model  was 
first  studied  as  early  as  in  19  53  by  H.  K.  Weiss  in  BRL  Re- 
port No.  879,  "Methods  for  Computing  the  Effectiveness  of 
Area  Weapons."   The  formula  for  the  salvo  kill  probability 
(2.7)  given  in  the  next  section  is  attributed  to  him.   Later, 
Breaux  [Ref.  7]  found  that  this  formula  was  not  suitable  for 
computation  when  n  is  large,  and  gave  another  method  of  com- 
putation.  However,  it  seems  that  there  is  no  published  re- 
search investigating  this  model  in  full  depth;  the  studies 
usually  come  to  an  end  when  expressions  for  the  kill  proba- 
bility are  derived.   In  this  chapter,  we  investigate  first 
the  property  of  the  kill  probability  as  a  function  of  param- 
eters involved  in  the  model.   Methods  for  computing  the  kill 
probability  and  bounds  for  the  kill  probability  are  also 
dealt  with  in  the  later  sections. 
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A.   SALVO  KILL  PROBABILITY  AS  A  FUNCTION  OF  n,  a  ,  a  ,  AND  a 

The  salvo  kill  probability  is  derived  in  two  steps  as 
is  stated  in  Chapter  I.   First  the  conditional  kill  proba- 
bility by  the  ith  weapon  is  calculated  given  that  the  random 
bias  (U,V)  takes  a  value  (u,v)  . 


CO       CO 


p(u,v)   =   /   /   D(u-x,v-y) f_(x,y)dxdy  .  (1.1) 


—  OO     — OO 


The  salvo  kill  probability  is  then  obtained  by  averaging 
l-{l-p(u/v)}   over  the  distribution  of  (U,V) ,  as 


CO       CO 

p   =   1  -  /    /   U-p(u,v)  }nf1(u/v)dudv  .  (1.2) 

—CO     —CO 

Using  the  assumptions  (2.2)  and  (2.3),  we  easily  calcu- 
late the  conditional  kill  probability  (1.1) 


/         >2W         ,2         2      2 
(u-x)    -i-(v-y)      _  x   +y 

,  0°  oo  2  „     2 

p(u,v)       =      -ir  /         /         e  2a  2°2    dxdy 

2  tt  a  2      -°°      -°° 


and  get 


2      2 

u  +v 

2  2(a2+oh 

p(u,v)       =      -£ j  e  z       .  (2.4) 

a     +a2 


17 


Substituting  (2.1)  and  (2.4)  into  (1.2),  the  salvo  kill 
probability  is 


2   2        2   2 

u  +v       u  +v 


2.2,       .2 


-,    °°   °°        2      2  (a  +a  )  n    2a, 
P   =   1    -^2"  /    /   tl  -  -^ J  e  }  e       dudv 


2ira-   -°°   -oo     a  +a~ 


Conversion  from  Cartesian  coordinates  into  polar  coordinates 
gives 


2  2 

r  r 


00       2      2(ct  +an)  n    2a. 


P   =   1  -  -\   /  {1  -  A j  e        2  }   e    ~1    r  dr  . 

a,   0     a  +a 


From  this  we  get 


P   =   1  -  /  (1  -  A  e"t/p)ne"tdt  (2.5) 

0      p 


where 


2   2  2    2    2 

X   =   a  /a  ,      p  =   (a  +o2)/o       .  (2.6) 


There  are  4  parameters  in  our  model,  but  the  salvo  kill 
probability  is  determined  by  three,  n,  a, /a,  and  cr2/a,  or 
n,  A,  and  p  as  is  seen  in  (2.5)  . 

If  the  integrand  in  (2.5)  is  expanded  in  a  binomial 
series  and  integrated  term  by  term,  we  obtain  the  Weiss 
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formula  for  the  salvo  kill  probability. 


P   =    l     (-U^MxV  4-  .  (2.7) 

-;  _i        J   p   J  tp 


P  is  obviously  an  increasing  function  of  n,  and  is  de- 

2   2 
creasing  in  a  /a  .   The  latter  is  easily  verified  by  in- 

2   2 
specting  (2.5),  where  A/p  is  constant  for  fixed  o?/a.    .   The 

2   2 

a^/a  -dependence  is,  however,  not  obvious.   In  (2.5)  ,  A  is 

2   2 
constant  for  a  fixed  value  of  a, /a  ,  and  there  are  two  factors 

which  involve  p  in  the  integrand,  with  p  tending  to  increase 

P  in  one  and  to  decrease  P  in  the  other. 

2   2 
Table  1  and  Table  2  are  given  here  to  show  the  a, /a  -  and 

2   2 
o~/a   -dependence  of  the  salvo  kill  probability,  where  n  is 

kept  constant,  n  =  2  for  Table  1  and  n  =  16  for  Table  2.   For 

2   2  2   2 

fixed  n  and  o./a  ,  P  is  decreasing  in  a   /a      as  mentioned 

2   2 
above.   However,  the  a./a  -dependence  is  a  little  different. 

2   2 
It  is  observed  that,  generally,  P  is  decreasing  in  a^/a      for 

2   2  2   2 

small  values  of  a, /a  ,  but  for  a  larger  a, /a  ,  P  increases 

2   2 
when  Op/a   is  increased,  reaches  a  maximum  value,  and  then 

decreases.   It  is  also  noted  that  from  the  viewpoint  of  real 

2   2 
world  applications,  the  case  of  large  a, /a  with  large  n  is 

especially  interesting;  e.g.,  as  seen  in  Table  2,  the  kill 

2   2 
probability  for  n  =  16  and  a,  /a   =  8.0  is  improved  to  0.5881 

2   2 
when  a^/a   is  increased  to  a  suitable  value.   The  maximum  P 

2   2 
is  more  than  1.7  times  the  original  value  for  a2/a   =  0. 
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Table  1.   Salvo  Kill  Probability 

n  =  2 


2  ,    2 


~2,  2 
a2/a 

0.00 

0.50 

1.00 

2.00 

4.00 

8.00 

0.0 

1.0000 

0.8333 

0.6667 

0.4667 

0.2889 

0.1634 

0.1 

0.9917 

0.8171 

0.6591 

0.4669 

0.2923 

0.1666 

0.2 

0.9722 

0.7977 

0.6487 

0.4647 

0.2940 

0.1689 

0.3 

0.9467 

0.7767 

0.6365 

0.4609 

0.2946 

0.1706 

0.4 

0.9184 

0.7550 

0.6232 

0.4560 

0.2944 

0.1717 

0.5 

0.8889 

0.7333 

0.6095 

0.4502 

0.2935 

0.1724 

0.6 

0.8594 

0.7120 

0.5956 

0.4439 

0.2920 

0.1728 

0.7 

0.8304 

0.6912 

0.5818 

0.4373 

0.2902 

0.1730 

0.8 

0.8025 

0.6712 

0.5681 

0.4305 

0.2881 

0.1729 

0.9 

0.7756 

0.6518 

0.5547 

0.4236 

0.2858 

0.1726 

1.0 

0.7500 

0.6333 

0.5417 

0.4167 

0.2833 

0.1722 
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Table  2    Salvo  Kill  Probability 
n  =  16 


2  ,    2 
a,  /a 


2/  2 

a2/a 

0.00 

0.50 

1.00 

2.00 

4.00 

8.00 

0.00 

1 .0000 

0.9935 

0.9412 

0.7835 

0.5511 

0.3370 

0.25 

1.0000 

0.9956 

0.9570 

0.8234 

0.6013 

0.3784 

0.50 

1.0000 

0.9965 

0.9659 

0.8504 

0.6399 

0.4128 

0.75 

1.0000 

0.9968 

0.9710 

0.8692 

0.6700 

0.4417 

1.00 

1.0000 

0.9967 

0.9739 

0.8824 

0.6938 

0.4662 

1.25 

0.9999 

0.9962 

0.9752 

0.8915 

0.7126 

0.4871 

1.50 

0.9997 

0.9954 

0.9753 

0.3978 

0.7276 

0.5049 

1.75 

0.9993 

0.9941 

0.9745 

0.9017 

0.7394 

0.5201 

2.00 

0.9985 

0.9922 

0.9728 

0.9038 

0.7485 

0.5330 

2.25 

0.9972 

0.9899 

0.9704 

0.9045 

0.7555 

0.5440 

2.50 

0.9954 

0.9869 

0.9672 

0.9038 

0.7606 

0.5534 

2.75 

0.9930 

0.9833 

0.9635 

0.9021 

0.7642 

0.5612 

3.00 

0.9900 

0.9792 

0.9591 

0.3996 

0.7665 

0.5678 

3.25 

0.9863 

0.9745 

0.9541 

0.8962 

0.7676 

0.5732 

3.50 

0.9821 

0.9692 

0.9486 

0.8921 

0.7677 

0.5776 

3.75 

0.9772 

0.9635 

0.9427 

0.8875 

0.7670 

0.5810 

4.00 

0.9719 

0.9573 

0.9364 

0.8824 

0.7655 

0.5837 

4.25 

0.9660 

0.9507 

0.9298 

0.8769 

0.7634 

0.5857 

4.50 

0.9597 

0.9438 

0.9228 

0.8710 

0.7607 

0.5870 

4.75 

0.9530 

0.9366 

0.9155 

0.8648 

0.7576 

0.5878 

5.00 

0.9459 

0.9291 

0.9081 

0.8583 

0.7541 

0.5881 

5.25 

0.9386 

0.9214 

0.9004 

0.8516 

0.7502 

0.5880 

5.50 

0.9309 

0.9135 

0.8926 

0.8448 

0.7460 

0.5874 

5.75 

0.9231 

0.9055 

0.8847 

0.8378 

0.7416 

0.5865 

6.00 

0.9151 

0.8973 

0.8767 

0.8307 

0.7369 

0.5853 
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Experience  shows  that  there  is  one  and  only  one  local  maximum 
of  P  as  a  function  of  a  /a   for  fixed  n  and  a, /a  ,  but  we 
have  not  succeeded  in  obtaining  a  proof. 

The  salvo  kill  probability  can  be  increased  by  increasing 

2   2  2   2 

o2/a   from  zero  for  given  n  and  a, /a   if  the  condition  below 

is  satisfied. 


Proposition  2.1 

The  salvo  kill  probability  P  for  fixed  n,  a   and  a  is 


not  maximal  at  o.  =  0  if  and  only  if 


+  oTT  +  ••  •  +  ~rr  >   1  f  (2-8) 


1+A    2  +  A  '  • • *  '  n+A 


where 


2/  2 

X      =      a   /a. 


Proof 

We  will  show  that 


dP/dp|      >   0  (2.9) 

p  =  A 


is  equivalent  to  (2.8).   Using  (2.7),  (2.9)  is  written  as 


n  /"(l-e-^V-V^^dt  <  n  /  (l-e-t/A)  ^V^^t/A  dt 
0  0 

(2.10) 
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Now  define 


I(x)   =   n  /  d.e-t/A^-l^t/A-xt  dt  _ 
0 


Then  the  inequality  (2.10)  is  rewritten  as 

1(1)   <   -I'  (1)/A  ,  (2.11) 

where  I1 (1)  means  the  value  of  the  derivative  dl/dx  at  x  =  1 
Let  s  =  exp(-t/A)  in  the  integral,  then  I (x)  is 


1  , 

I(x)   =   nA  /   sAX(l-s)n_1  ds 


=   A    n    n"1 


Ax+1  Ax+2      Ax+n 


Therefore, 


I'(x)         y    A 
I(x)         ^   Ax+j 


The  inequality  (2.8)  is  readily  obtained  from  this 
equation  with  x  =  1  and  (2.11) .  □ 

Proposition  2.2 

A  sufficient  condition  that  the  salvo  kill  probability 
P  for  fixed  n,  a, ,  and  a  is  not  maximal  at  o.  =  0  is 
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a2/a^   ^   n/(e-l)  -  1  .  (2.12) 


Proof 

Since  l/(x+A)  is  decreasing  in  x, 


n    i        n    J+1  ^        n+1   i  n,ul 

;Z!  5+A      •  i   •    X+A      i    X+A  1+A 


If  the  right-hand-side  is  equal  to  or  greater  than  unity, 
i.e.,  l+n/(l+A)  =  e,  then  the  condition  (2.8)  in  Proposition 

2.1  is  satisfied.   The  relation  (2.12)  is  immediate  from  the 

above-mentioned  inequality.  G 

Proposition  2.3 

Assume  that  (2.8)  holds.   A  necessary  condition  that  the 

salvo  kill  probability  be  maximum  at  some  point  a~  >  0  for 
fixed  n,  a,  and  a  is 


I    (-1)  j"1(^)j(^)j"1  p+^'12      -   °>  (2-13) 


j=l        J    w  (p+jl 


where 


A   =   a  /a,     ,  p   =   (a  +  a2)/a^  . 


Proof 

Differentiate  (2.7)  with  respect  to  p,  and  equat  it  to 
0.   The  equation  (2.13)  is  immediate.  D 
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Example.   Case  n  =  2 

Equation  (2.13)  is  generally  too  complicated  to  be  solved 
for  cr_/a  in  closed  form.   Here  we  solve  it  for  the  simplest 
case,  n  =  2. 

The  inequality  (2.10)  with  n  =  2  is 


1  *4r  >   1  , 


1+A    2+A 


and  after  a  bit  of  algebra,  this  reduces  to 


1/A   =  q\/o}      >      2/C/5  -  1)  .  (2.14) 


2   2 

Now,  assume  that  a, /a   is  large  and  satisfies  (2.14).   From 

(2.13) 


_£ -  2  —   P  +  1 


=   0 


(p+D2      P  (p+2)2 
Therefore, 


2/2       ,,    ,       (p+1)  3 
a9/a   =   p/A  -  1   =  — ^ '-j   -  1  , 

p(p+2)Z 


r,2/    2             i  /y                   (p+D3 
a-, /a   =   1/A   =  —5 — j   ' 

p  (P+2P 


(2.15) 


2  ,  2 


This  is  a  parametric  expression  of  maximizing  a~/a      as  a 

2   2 
function  of  a-,  /a   the  curve  of  which  is  shown  in  Fig.  2.   The 

value  of  the  parameter  p  which  corresponds  to  the  point 
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o 
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{o2/a2   =    (/5+D/2,    a2/a2   =    0)    is    p   =    pQ,    where 


Pq      =      (/5   -    l)/2    . 


When  p  approaches  zero, 


o2/a2      -      1/4P2  ,        a2/a2   -   1/4 p  , 


and  therefore,  the  asymptotic  form  of  the  curve  is 


2   2      1 
a2/a   =  y  a, /a  .  (2.16) 


B.   COMPUTING  THE  SALVO  KILL  PROBABILITY 

In  the  preceding  section,  the  Weiss  formula  for  the  salvo 
kill  probability  was  given  in  (2.7).   According  to  Breaux 
[Ref.  7],  this  formula  is  not  suitable  for  calculation  for  a 
large  n  because  it  is  an  alternating  series.   Breaux  pointed 
out  that  the  salvo  kill  probability  could  be  expressed  in 
terms  of  an  Incomplete  Beta  Function,  and  that  for  large  n, 
Tang's  method  [Ref.  8]  would  be  advisable  for  calculating  an 
Incomplete  Beta  Function. 

First,  we  present  a  recurrence  formula  for  calculating  the 
salvo  kill  probability  which  is,  in  principle,  the  same  as 
that  of  Tang's  method.   Its  generalization  will  be  given  in 
a  later  chapter. 
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Proposition  2.4 

Let  us  denote  the  salvo  miss  probability  by  n  weapons  as 
Q(n)  =  1-P. 


00 


Q(n)   =   /  (l-ye~t/p)n  e"t  dt  .  (2.17) 

0 


Then  Q(n)  satisfies  the  following  recurrence  equations: 


Q(n)   =   HTp-  Qtn-1)  +  -E_(i-P)"  ,    n  =  1,2,...        (2.18) 


Q(0)   =   1  ,  (2.19) 


where 


2    2    2  2    2    2 

y   =   A/p   =   a  /  (a  +a2)  '    P   =   Ca  +0^)/^^ 


Proof 

From  (2.17) 


Q(n)   =   /  (l-ue"t/p)n  e_t  dt 


r  /i    "~t/p.n-l  -t,.    t    ,,    -t/p.n-^1  -t/p-t,. 
J  (1-ye  '  w)         e   dt  -  yj  (1-ye    K)    e      dt 

0  0 
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The  first  term  in  the  right  hand  side  is  Q(n-l).   In  the 
second  term,  we  integrate  it  by  parts,  and  have 

00  oo 

Q(n)   =   Q(n-l)  -  jMl-ye     )  e   |   -  £-  J     (1-ye   '  )  e 


dt 


=   Q(n-l)  +  fd-u)n  "  £  Q(n)  . 


The  equation  (2.18)  is  immediate.  G 

It  is  worthwhile  to  note  that  the  method  presented  in 
Proposition  2.4  gives  not  only  the  value  of  Q(n)  for  a 
specified  n  ,  but  also  gives  all  the  Q(n) '  s  up  to  n_ .   It 
is  an  important  advantage  of  this  method  over  the  Weiss 
formula  because  we  are  often  interested  in  the  salvo  kill 
probabilities  for  several  n's.   As  is  shown  in  Table  3,  the 
computational  time  for  Q(n)  by  our  recurrence  algorithm  is 
about  twice  the  time  required  for  computation  of  a  single 
Q(n)  using  the  Weiss  formula.   When  the  salvo  kill  probabili- 
ties for  more  than  a  single  n  are  needed,  therefore,  use  of 
the  recurrence  algorithm  given  in  Proposition  2.4  is 
recommended . 

We  have  another  efficient  algorithm  for  computing  the 
salvo  kill  probability,  which  is  good  even  for  large  n. 
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Table  3.   Comparison  of  the  Processing  Time  for 
Computing  the  Salvo  Kill  Probability 


Mean  time  (microsecond),   p  =  1,  and  \i   =  0.5. 


n 

10 
20 
30 


Weiss  formula 

92 
189 
300 


Recurrence  formula 

175 
383 
587 


Mean  time  (microsecond),   p  =  1,  and  n  =  10. 


0.1 

0.5 
0.9 


Weiss  formula 

90 

92 

100 


Recurrence  formula 

196 
175 
186 


Computation  is  on  the  IBM  30  33  installed  in  the  Naval 
Postgraduate  School.   Each  value  is  the  mean  of  4  trials. 
One  trial  is  1000  repetitions  of  calculation,  and  the 
elapsed  time  units  are  divided  by  10  00  and  multiplied  by 
2  6  (microseconds)  . 
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Proposition  2.5 

The  salvo  miss  probability  Q(n)  by  n  weapons 


Q(n)   =   /  (1  -ye~t/p)n  e_t  dt  (2.17) 

0 


satisfies  the  following  relationship 


n  ,   , 

Q(n)   =    I       (?)  (l-y)n"K  yK  q(k)  ,  (2.20) 

k=0   K 


where  q(k)  is  the  salvo  miss  probability  by  k  weapons  with 
V  =  1, 

00 

q(k)   =   /   (1  -e~t/p)k  e_t  dt  . 
0 


Proof 


Q(n)   =   /  Ul-y)  +  y(l-e~t/p)}n  e_t  dt 
0 


r   ,n.  ,,   v  n-k   k  t       ,,     -t/p.k   -t  ,. 
I      (,)  (1-y)     \i     J       (1-ye  /H)   e   dt 

k=0   K  0 


D 
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Corollary  2.6 

The  salvo  miss  probability  (2.17)  can  be  calculated  by 
the  formula 


n 
Q(n)   =    I      A.  ,  (2.21) 

k=0   K 


where 


Ak/Ak-1   =   u(n-k+l)/(l-y)  (k+p)  ,   k=l,2,... 

(2.22) 
AQ   =   (l-y)n  . 

Proof 

As  can  be  seen  from  (2.18)  and  (2.19)  with  y  =  1,  q(0) 
is  unity,  and  q(k)/q(k-l)  is  k/(k+p).   Utilizing  the  fact 
that  (v)/(j__t_)  =  (n-k+l)/k,  the  corollary  can  be  obtained  by 
substitution.  D 

Note  that  all  the  terms  in  the  expression  (2.24)  are 
positive,  so  there  is  no  problem  arising  from  cancellation 
of  terms  with  alternating  signs  as  was  observed  in  the  Weiss 
formula  (2.7)  . 

In  checking  the  computation,  the  following  P's  for  special 
cases  would  be  of  use. 


(1)   The  case  a,  =0  corresponds  to  p  ■*■   «,  and  we  have 
P   =   1  -  (1  -y)n  . 
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(2)   The  case  a      =  0  corresponds  to  y  =  1 ,  p  =  A ,  and  so 
from  the  recurrence  formula,  we  have 


p   -  i       n(n-l)  •  •  •  1 

(n+A)  (n-l+A)  •••  (1+A)  *  \n<k) 


2    2     2 
(3)   The  case  a  +a2  =  a,.   The  formula  (2.5)  with  p  =  1 


is  easily  integrated  to  give 


P   =   1  -  jjij-U  -  d-A)n+1}  ,  (2.25) 


where 


2  ,  .    2  2,  .  2,2 

p   =   a  /(a   +  a2)  ,    A   =   a  /a,  . 


C.   BOUNDS  TO  THE  SALVO  KILL  PROBABILITY 

In  Chapter  I  it  was  pointed  out  that  an  upper  bound  of 
the  salvo  kill  probability  is  given  by  the  kill  probability 
of  n  repeated  independent  shots. 

P   <   P 


where 


00       00 


=      1   -    [/        /      {1  -p(u,v)  }f1(u/v)dudv]n    . 


—oo    — oo 


Under  the  assumptions  of  the  present  chapter,  we  have 
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Proposition  2.7 

An  upper  bound  to  the  salvo  kill  probability  is  given  by 


P   =   1  -  (1  -  j^)n    ,  (2.26) 


where 


^       2/~2  /  2  J  2W  2 

A   =   a  /Oj^  ,    p   =   (a  +  02)/^ 


As  to  a  lower  bound,  we  have 


Proposition  2.8   (Merritt  and  King) 

A  lower  bound  to  the  salvo  kill  probability  is  given  by 


P   =   1  -  p(ny)~p  Tnii(p)  ,  (2.27) 


where  r  (p)  is  the  Incomplete  Gamma  Function, 


r  (p)   =   /   e  r  tp  L   dt  , 
x        0 


and 


y   =   a  /(a  +  a,)  ,    p   =   (a  +o2)/o->     . 


Table  4  and  Table  5  show  how  the  true  kill  probabilities 
are  bracketed  by  these  two  bounds.   The  triplets  of  entries 
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Table  4    Salvo  Kill  Probability  and  Its  Bounds 

n  =  2 


2  . 

a1/a 

2 

2  ,    2 
a2/a 

0.00 

0.50 

1.00 

2.00 

4.00 

3.00 

0.3889 

0.7500 

0.6400 

0.4898 

0.3306 

0.1994 

0.5 

0.8889 

0.7333 

0.6095 

0.4502 

0.2935 

0.1724 

0.7364 

0.6187 

0.5216 

0.3911 

0.2580 

0.1527 

0.7500 

0.6400 

0.5556 

0.4375 

0.3056 

0.1900 

1.0 

0.7500 

0.6333 

0.5417 

0.4167 

0.2833 

0.1722 

0.6321 

0.5443 

0.4715 

0.3679 

0.2532 

0.1552 

0.6400 

0.5556 

0.4898 

0.3951 

0.2840 

0.1814 

1.5 

0.6400 

0.5524 

0.4825 

0.3829 

0.2696 

0.1689 

0.5507 

0.4832 

0.4267 

0.3429 

0.2442 

0.1542 

0.5556 

0.4898 

0.4375 

0.3600 

0.2653 

0.1736 

2.0 

0.5556 

0.4881 

0.4333 

0.3524 

0.2554 

0.1643 

0.4866 

0.4333 

0.3882 

0.3192 

0.2338 

0.1516 

0.4898 

0.4375 

0.3951 

0.3306 

0.2489 

0.1664 

2.5 

0.4898 

0.4365 

0.3925 

0.3255 

0.2418 

0.1593 

0.4353 

0.3922 

0.3554 

0.2977 

0.2233 

0.1482 

0.4375 

0.3951 

0.3600 

0.3056 

0.2344 

0.1597 

3.0 

0.4375 

0.3944 

0.3583 

0.3021 

0.2292 

0.1542 

0.3935 

0.3580 

0.3274 

0.2784 

0.2131 

0.1444 

0.3951 

0.3600 

0.3306 

0.2840 

0.2215 

0.1536 

3.5 

0.3951 

0.3596 

0.3294 

0.2815 

0.2175 

0.1492 

0.3588 

0.3291 

0.3032 

0.2611 

0.2034 

0.1404 

0.3600 

0.3306 

0.3056 

0.2653 

0.2099 

0.1479 

4.0 

0.3600 

0.3303 

0.3048 

0.2635 

0.2068 

0.1443 

0.3297 

0.3045 

0.2823 

0.2457 

0.1944 

0.1365 

Triplets  are,  from  the  top,  P,  the  upper  bound, 
P,  the  true  value,  and  P,  the  lower  bound. 
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Table  5   Salvo  Kill  Probability  and  Its  Bounds 

n  =  16 


2  ,  2 
ax/a 


2.    2 

a2/a 


1.0 


2.0 


3.0 


4.0 


5.0 


6.0 


7.0 


8.0 


0.00 

1.0000 
1 .0000 
0.9997 

0.9985 
0.9885 
0.9952 

0.9900 
0.9900 
0.9817 

0.9719 
0.9719 
0.9592 

0.9459 
0.9459 
0.9305 

0.9151 
0.9151 
0.8983 

0.8819 
0.8819 
0.8647 

0.8481 
0.8481 
0.8310 


0.50 

0.9997 
0.9967 
0.9944 

0.9954 
0.9922 
0.9862 

0.9821 
0.9792 
0.9685 

0.9597 
0.9573 
0.9432 

0.9309 
0.9291 
0.9130 

0.8987 
0.8973 
0.8802 

0.8650 
0.8639 
0.8468 

0.8313 
0.8305 
0.8136 


1.00 

0.9985 
0.9739 
0.9688 

0.9900 
0.9728 
0.9644 

0.9719 
0.9591 
0.9469 

0.9459 
0.9364 
0.9215 

0.9151 
0.9081 
0.8917 

0.8819 
0.8767 
0.8597 

0.8481 
0.8441 
0.8273 

0.8147 
0.8117 
0.7953 


2.00 

0.9900 
0.8824 
0.8750 

0.9719 
0.9038 
0.8935 

0.9459 
0.8996 
0.8864 

0.9151 
0.8824 
0.8674 

0.8819 
0.8583 
0.8424 

0.8481 
0.8307 
0.8145 

0.8147 
0.8016 
0.7857 

0.7824 
0.7724 
0.7570 


4.00 

0.9459 
0.6938 
0.6867 

0.9151 
0.7435 
3.7388 

0.8819 
0.7665 
0.7546 

0.8481 
0.7655 
0.7522 

0.8147 

0.7541 
0.7402 

0.7824 
0.7369 
0.7229 

0.7515 
0.7167 
0.7030 

0.7222 
0.6952 
0.6819 


3.00 

0.8147 
0.4662 
0.4611 

0.7824 
0.5330 

0.5253 

0.7515 
0.5678 
0.5590 

0.7222 
0.5837 
0.5739 

0.6945 
0.5881 
0.5778 

0.6684 
0.5853 
0.5749 

0.6439 
0.5780 
0.5677 

0.6209 
0.5680 
0.5579 


Triplets  are,  from  the  top,  P,  the  upper  bound, 
P,  the  true  value,  and   P,  the  lower  bound. 
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in  the  table  are,  from  the  above,  the  upper  bound  P,  true 

value  P,  and  P,  the  lower  bound.   It  is  observed  that 

—  '   2 

P  -  P  is  not  small  for  smaller  values  of  o-/a    ,  but  gets  smaller 

2   2  2   2  — 

when  a   /a      is  increased.   If  a, /a   is  small,  P  and  P  are 

close.   On  the  other  hand,  P  is  fairly  close  to  P  over  wide 

2   2       2   2 
ranges  of  a, /a   and  o~/a    .   Figures  3  and  4  illustrate  this. 

It  seems  reasonable  that  Merritt  and  King  used  P  as  an 

approximate  formula  for  P  in  their  study  on  the  optimal 

ballistic  dispersion,  because  computation  was  really  a 

problem  at  that  time.   Nowadays,  however,  direct  computation 

of  P  is  much  more  natural. 
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1.00 


0.98  - 


0.96  - 


0.94 


n=  20 

— 

V         P 

**/a2  =  1.0 

—  / 
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1 

1 

1              \ 
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2.0  3.0 
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4.0 


5.0 


Figure    3^. 
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III.   OPTIMIZATION  OF  SALVO  FIRING 

In  Chapter  II  we  have  learned  the  following:   The  salvo 
kill  probability  is  a  decreasing  function  of  o.  if  other 
parameters  are  kept  constant.   Therefore,  to  obtain  higher 
salvo  kill  probability,  it  is  necessary  to  make  a,  as  small 
as  possible.   The  salvo  kill  probability  as  a  function  of 
a2  is,  however,  complicated.   If  a,    is  sufficiently  small, 
the  salvo  kill  probability  is  monotone  decreasing  in  a~, 
but  if  a,  is  not  so  small  and  satisfies  the  condition  (2.8) , 
then  the  salvo  kill  probability  is  not  maximal  at  a.  =  0, 
but  takes  a  maximum  value  at  some  a„  >  0.   In  other  words,  the 
smallest  ballistic  dispersion  is  not  always  most  desirable. 
Deliberate  increase  of  the  ballistic  dispersion  is  encoun- 
tered in  shotguns  and  in  related  military  weapons.   In  this 
chapter  we  present  an  approximate  formula  for  the  optimal 
ballistic  dispersion.   We  will  first  briefly  sketch  Walsh's 
theory  [Ref.  9]  since  our  approximate  formula  is  based  on 
it.   In  Section  B  our  approximate  formula  is  presented. 
Its  accuracy  will  be  studied  in  Chapter  IV. 

A.   WALSH  THEORY 

The  Walsh  model  is  based  on  assumptions  similar  to  those 
of  Chapter  I.   A  salvo  of  n  weapons  is  fired  at  a  target. 
The  impact  point  error  is  composed  of  two  parts,  a  random 
bias  which  is  common  to  all  n  weapons,  and  a  round-to-round 


40 


error.   A  Cartesian  coordinate  system  is  chosen  so  that  its 
origin  coincides  with  the  center  of  impact  of  n  weapons. 
The  coordinates  of  the  random  position  of  the  target  are 
denoted  as  U  and  V  with  respect  to  this  system.   The  joint 
probability  density  of  U  and  V  is  given  by  f , (u,v) .   The 
impact  points  of  weapons  (X.,Y.)  are  independent  and  iden- 
tically distributed  random  variables  with  joint  density 
f2(x,y) .   It  is  assumed  that  impact  of  a  weapon  at  (x,y) , 
given  the  target  is  at  (u,v) ,  destroys  the  target  with 
probability  D(u-x,v-y)  independently  of  other  weapons. 
The  lethal  area  of  a  weapon  is  denoted  as  A: 

00       CO 

j         f        D(x,y)&xdy   =  A  . 

—00     —CO 

The  salvo  kill  probability  is  given  by  (1.2), 


00       00 

P   »   1  -  /    /{I-  p(u,v) }nf1(u/v)dudv  ,        (1.2) 

— OO     —CO 

where  p(u,v)  is  the  conditional  destruction  probability  by 
a  single  weapon  given  that  the  target  is  at  (u,v) . 


OO       OO 

p(u,v)   =   /    /    D(u-x,v-y) f2(x,y) dxdy  .         (1.1) 

—CO      — OO 

Note  that  p(u,v)  can  be  integrated  to  give  A. 
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00       00  00       00 

/    /    p(u,v)dudv   =   A  /    /    f2(x,y)dxdy   =   A  .   (3.1) 

— OO     —CO  —00     —00 


The  problem  is  to  derive  the  optimal  ballistic  disper- 

2 
sion,  a„.   Instead  of  solving  the  original  problem,  Walsh 

solved  a  revised  problem  that  was  different  from  the  original 

2 
one  in  that  the  variable  is  not  a~  itself  but  rather  the 

function  p(u,v) . 
Proposition  3.1  (Walsh) 


00  OO 

max  P      =      1    -   /         /         {l-p(u,v)}    f,  (u,v)dudv  (3.2) 

pt-0 


— OO  —CO 


subject   to 


00  00 


/         I         p(u,v)  du  dv       =      A    , 


—CO  —CO 


p(u,v)       =      0 


has  a  solution 


where 


(3.3) 


1  -  {fYiW}n"1     (u'v)  £  E 

13.41 

(u,v)     6    E 
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E   =   { (u,v) |f1(u/v)  1   C)  .  (3.5) 


C  is  positive  and  is  given  by 


/  /  p*(u,v)  dudv   =   A  .  (3.6) 

E 


The  corresponding  salvo  kill  probability  is 


P*   =   CA  +  /  /  {f,  (u,v)  -  C}  dudv  .  (3.7) 

E 


In  case  the  random  target  bias  is  given  by  a  normal 
distribution,  we  also  have 

Corollary  3.2 


00       00 

max  P   =   l-/    /   {1  -  p(u,v)  }nf,  (u,v)  dudv       (3.8) 
p(") 


—CO     —CO 


subject  to 


CO       CO 


(        j        p(u,v)  dudv   =  A  , 


—CO     —CO 


(3.9) 


p(u,v)   >   0  , 


where 


2     2 
u    v 


0  2   0  2 

f,  (u,v)   =   ^ e     U    V  (3.10) 

1  27Tau°v 
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has  the  solution 


2    2 


1  -  exp{-4>  +  2  (n-.i)  (^2"  *  ^2^   (u'v)  e  E 

U   V  (3.11) 

0  (u,v)  e  E 


where 


E   =   {(u,v)|u2/a2  +  v2/a2  <  2(n-l)<f>} 

1       U  V  "" 


$  is  positive  and  is  determined  by 


e"*  -  1  +  <j>  -  A/2ir(n-l)a  a    =   0  .  (3.12) 

u  v 


The  maximum  value  of  the  salvo  kill  probability  is  given  by 

P*   =   1  -  e"(n""1)(|){l  +  (n-1)  (.1  -e"*)  }  .  (3.13) 

Proof 

The  equation  (3.11)  is  readily  obtained  from  (3.4)  in 
Proposition  3.1  with  f, (u,v)  given  by  (3.10)  and  C, 


C   =   ,-J: e-(n-D(J)  (3el4) 

2tt  a  a 
u  v 


Substituting  (3.11)  in  (3.6)  and  changing  variables  into 
(r,6)  by 
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u      =      /2(n-l)   aur  cos    0  ,       v      =      /2(n-l)   a     r  sin    6  (3.15) 


we   have 


A      =      /    /    P*(u,v)  dudv 
E 


=      2ir(n-l)o     a       /         (1  -  e-<p    r    )  2r   dr 

0 


=      27T(n-l)  auav  (<j>  -1  +  e"*)     . 


Substituting    (3.10)     and    (3.14)     in    (3.7)    with    (3.15), 

*  A  -(n-l)<j>   ,        2tt       ,       .  ,  (    *         ,      ..     2  ,       ,,  , 

P*      =      0^   -      e  +t (n-l)a   a      /       ,    -(n-l)r        -   n-li,    a 

2ira   a  2tto    a  u   v  ;.     (e  -e  )  2rdr 

=      1    -    e"(n"1)(}>{l    +    (n-l)cj)    -   A/2ira     a     }     . 

u    v 

This  formula  together  with  (3.12)  gives  (3.13).       □ 

For  simplicity,  let  us  call  the  function  p*(u,v)  the 
(Walsh)  optimal  p  function,  and  the  associated  salvo  kill 
probability  P*  the  Walsh  optimal  P.   It  should  be  noted  that 
the  constraint  in  the  revised  problem  is  weaker  than  that  of 
the  original  problem,  and  so  the  Walsh  optimal  P  is  an  upper 
bound  to  the  maximum  salvo  kill  probability  of  the  original 
problem.   It  should  also  be  noted  that  the  Walsh  optimal  p 
function  is  not  feasible  in  the  original  problem  (being  non 
zero  on  only  a  finite  set)  if  dispersion  is  present. 
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B.   APPROXIMATE  FORMULA  FOR  THE  OPTIMAL  BALLISTIC  DISPERSION 

Walsh  developed  the  theory  stated  above,  but  he  did  not 
extend  his  theory  beyond  the  results  (3.4)-(3.7).   The  impor- 
tant part  left  by  Walsh  and  subsequent  workers  is  to  bridge 
the  gap  between  the  revised  problem  and  the  original  prob- 
lem.  Walsh  got  the  optimal  p  function  as  given  in  (3.4) , 
whereas  we  need  the  optimal  density  function  f2(x,y) .   The 
material  in  this  section  bridges  that  gap. 

Let  us  define  three  double  Fourier  transforms  as  follows: 


00    °°    i  (go,  u+go^v) 
p*(u)^,u>2)   =   /    /    e  p*(u,v)dudv, 

—CO    — oo 


00    °°    i  (oj,  x+Go^y) 
D  (a)., u>2)   =   /    /    e  D(x,y)dxdy,  (3.16) 

— oo    — oo 


00   °°    i(Go,x+co^y) 
f2(Go1,Go2)   =   /    /    e  f2(x,y)dxdy 

— oo    — oo 


Then  from  (1.1)  we  have 

P*(go,,go2)   =   D(go,  /0j2)  f  2  (go,  /Go2)  .  (3.17) 

In  the  Walsh  theory,  D(x,y)  is  a  given  function,  and 

the  Walsh  optimal  p  function  is  given  by  Proposition  3.1. 

Therefore,  both  D(go1,go2)  and  p*(go..,oj  )  in  (3.17)  are  now 
known  functions.   Thus  the  Fourier  transform  of  the  "optimal" 
f2(x'Y)  can  ke  calculated  by 
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P*  (oj^oO 


f2(a)1/a)2)   = = —  ,  (3.18) 


D(w,  ,(0 


and  we  will  get  f 2 (x,y)  as  its  inverse  Fourier  transform. 
This  might  have  been  the  end  of  Walsh  theory,  but  there  are 
difficulties.   First,  we  do  not  have  any  guarantee  that 
f2(x,y)  is  a  density  function,  because  Walsh  solved  the 
restricted  optimization  problem  subject  to  p(u,v)  >  0  but  not 
f„(x,y)  >  0.   Further,  even  if  this  problem  were  solved,  there 
would  still  remain  the  problem  of  feasibility  of  the  f„  func- 
tion.  Walsh  suggests  in  his  paper  that  it  would  be  possible 
to  obtain  a  real-world  dispersion  function  f2(x,y)  as  the 
theory  requires,  at  least  approximately,  by  careful  design 
of  the  ammunition,  but  this  may  be  asking  too  much  of  the 

manufacturing  process.   We  assume  below  that  the  function 

2 
f2(x,y)  is  normal,  with  only  the  variance  a2  subject  to 

manufacturing  control. 

The  problem,  then,  is  to  determine  or  approximate  the 

optimal  ballistic  dispersion  a2,  given  the  damage  function 

D(x,y)  and  the  function  p*(x,y)  from  (say)  Proposition  3.1. 

In  the  following,  both  D(x,y)  and  p*(x,y)  are  assumed  to  be 

even  functions  of  x  and  y  respectively,  and  to  have  "moments" 

of  all  degrees.   The  lethal  area  is  denoted  as  A. 

oo     oo  oo     oo 

/    /    xD(x,y)dxdy  =  /    /    yD(x,y)dxdy  =  •••  =  0  ,   (3.19) 

—  OO    —00  — oo    — oo 
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00  CO  00  oo 


/         /         xp*(x,y)dxdy  =    /         /        yp*(x,y)dxdy   =    •••    =    0  , 


— oo         — oo  — oo  — oo 


OO  00 


/         /         D(x/y)   dxdy      =      A    ,  (3.20) 


— OO  —00 


°°  °°  2k 

/         /        x      D(x,y)dxdy      <      •*>,       k  =   1,2,    .. 


— OO  — OO 


°°         °°  2k 

/         /        x      p*(x,y)   dxdy      <      «,    k  =    1,2,     ... 

—00  — OO 


etc. 

For  the  sake  of  generality,  we  actually  consider  the  non- 
symmetric  case  where  the  optimal  dispersion  may  be  differ- 
ent in  the  x  and  y  directions.   Let  the  two  dispersions  be 

2       2 
a      and  a  .   Then 
x      y 


1 

2 

X 

o    2 

2a 

e        x 

2 

.  Y 

o    2 
2a 

y 

27Tax  ay 

/ 

f9(x,y)   =   -   *  -   e   A  *    ,  (3.21) 

2    J       2Tra_.  a 


with  the  double  Fourier  transform 


2  2    2  2 

a  to,   a  u. 
x  1    y  2 

2   "   2 
f2(ai1,a)2)   =   e  .  C3.22} 
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Expanding  the  exponential  function  in  the  double  Fourier 
transform  of  D(x,y),  and  integrating  term  by  term,  we  have 


OO  00 


i  ((jo,  x+co2y) 
D(o>,  ,0)-)    =    /         /         e  D(x,y)dxdy 


— OO  — OO 


00  ,  -i  oo  oo 

I      YT  {        j         (oo1x+oj2y)  -"d  (x,y)   dx  dy 

"1=0      -1  —oo  — oo 


2 

=      Ml-yj/         /         x2D(x,y)dxdy 


—  OO  — oo 


2 

s2 
2    A 


2     1  2 

-   -5-  t  j        !        y   D(x,y)   dxdy   +    ...] 


—00      —00 


where   the  third   line    is   obtained  because   D   is   even.      There- 

2  2 

fore,    for    sufficiently   small   values   of   u),    and  co„, 


D(u>,,w2)       =      Dn^C0i'(;02^  (3.23) 


where 


2 
D0(co1,(Ja2)       =      A  exp  [ — j-  ^  J         f        x   D(x,y)   dx  dy 


—00      —00 


2 
2    A 


±.  /        /        yZD(x,y)   dxdyj     .  (3.24) 


—00      —00 


Similarly, 
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p*(a)1,a>2)       =      p0(co1,u2)  (3.25) 


2  2 

for   small  u>,    and  u>    ,    where 

2 

Wl    1      °°        °°         2 
Pgdo^u^)    =  A  exp  [ j-  -  /         /         x   p*(x,y)   dx  dy 

—  .jo  —00 


2 

yi/         /         yp*(x,y)   dxdy]     .  (3.26) 

—00      —00 


2       2 
Now  we  determine  two  parameters  a   and  a   in  (3.21)  by 

the  relation 


p0(a,1,w2) 

f2(u>,,u>2)   = •  (3.27) 

D0((0iru2) 

Then  from  (3.22),  (3.24),  and  (3.26)  we  get  an  approximate 
formula  for  the  optimal  ballistic  dispersion: 


00   00 


al   =  [i  /    /    x2{p*(x,y)-D(x,y)}dxdy,  0]  + 


x    lA 

—00    —00 


(3.28) 


00    00 
al   ■  [7  /    /    y  (p*(x,y)-D(x,y)}  dxdy,  0] 

■*  — c 


•00    —00 


where 


if   a    0 


( 

[a,0]+   = 

(  0 


otherwise 
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Note  1 

In  this  proposition,  the  approximate  formula  for  the 
optimal  dispersion  is  derived  from  the  condition  that 

f-tco-j/U)-)  coincides  with  f  _  (to..  ,tO  for  inf  initesimally  small 

2       2 
to,  and  to-.   It  would  be  better  if  we  could  derive  the  approxi- 
mate formula  from  the  condition  that  f-(.co,  ,gO  is  nearly 

~*  2       2 

equal  to  f2(to,,to2)  over  a  wider  domain  of  to,  and  oo2  which 

include  (0,0)  .   In  the  next  chapter  we  will  give  an  idea  of 

improving  the  approximate  formula  taking  this  into  account. 

Note  2 

Generally,  approximation  (3.23)  with  (3.24)  is  accurate 
only  up  to  the  second  power  of  to,  and  co_.   If  the  damage 
function  is  of  Gaussian  type, 


2     2 


D(x,y)   =   e 


0  2   0  2 

2a    2a 
x     y 


then 


DqCw^oO   =   D((u1,a)2) 

is  an  exact  relationship. 

As  a  corollary,  we  have  the  following  approximate  formula 
for  the  model  with  normal  errors  and  Gaussian  damage  function, 

Consider  a  salvo  firing  of  n  weapons.   The  density  func- 
tions of  the  random  bias  and  the  round-to-round  ballistic 
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error   are,    respectively, 


2  2 

u  v 


o    2      o    2 

1  2a        2a 

fl(u'v)      ■     2i,q    a      e        "  V    •  (3-32) 

U      V 


and 


- 

2 

X 

2 
.    Y 

n       2 

1 

e 

2a 

X 

2a 

y 

2tt 

a 

X 

a 

y 

/ 

f2(x^>      =      2F5    a      e  '  (3-33) 


the  damage  function  is  given  by 

2        2 

_ (u-x)    (v-y) 

o  2      .  2 
2a      2a 

D(u-x,v-y)   =   e     u       v  (3.34) 


Then  approximate  formulae  for  the  optimal  ballistic  dispersion 

for  fixed  n,  a    ,    a    ,    a    ,    and  a   are 
u   v   x       y 


(ox/ax)2  =  [i(n-l)2(|)2(au/ax)3((t)v/a  )-(n-l)  (au/ax)2-l,0]  + 


C3.35) 


(a  /a  )2  =  [i(n-l)  2$2  (o   /a  )3(a  /a  )-(n-l)(a  /a  )2-l,0J  + 
Y      Y  2  v   y    u   x         v   y 


where  4>  is  determined  by  the  equation 


e"*  -  l  +  cj,  -  (a  /a  )  (a  /a  )/(n-l)   =   0  .        13.36) 
T      x   u    y   v 
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An  upper  bound  to  the  maximum  salvo  kill  probability  is  given 
by 


=  1  -  e-(n-1)(^{l  +  (n-l)U-e-*)}  . 


(3.37) 


The  solution  to  the  Walsh  problem  with  the  above  men- 
tioned assumptions  is  already  given  in  Corollary  3.2.  By 
use  of  (3.11),  we  have 


OO  00 


j^  j        f        u  p*(u,v)  du  dv 


— »oo  —00 


3  r2 

auav   -      /2(n-1)(Jl   1     2  M  -*+37n=IJ\       - 

— r —  2ttJ  7r     11  -  e  1  rdr 

A       0       ^ 
S  v{(n-l) V  -  (n-l)A/7ra  a  } 


and 


OO       OO 


^  /    /    x  D(x,y)  dxdy   =   a^  . 


—00     — OO 


Substituting  these  relations  into  (3.28),  we  get  the  approxi- 

2 
mate  formula  for  a   (3.35).   Similar  arguments  follow  for 

2 

a    . 

y 
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IV.   SIMPLE  SALVO  MODEL  REVISITED   • 

In  this  chapter  we  apply  the  general  approximate  formula 
for  the  optimal  ballistic  dispersion  obtained  in  the  last 
chapter  to  the  Simple  Salvo  Model  dealt  with  in  Chapter  II 
and  investigate  its  accuracy. 

The  assumptions  in  this  chapter  are  the  same  as  those  in 
Chapter  II.   We  consider  a  salvo  of  n  weapons  against  a 
target.   The  impact  point  error  is  composed  of  two  parts, 
the  random  bias  and  the  round-to-round  error.   The  position 

of  the  target  (U,V)  has  a  circular  normal  distribution  centered 

2 
at  (0,0)  with  variance  a, .   The  impact  points  of  the  ith 

weapon  (X,Y)  are  assumed  independent  and  identically  dis- 
tributed circular  normal  random  variables  centered  at  (0,0) 

2 
with  variance  a2.   Finally,  we  assume  the  circular  Gaussian 

damage  function  with  parameter    : 


2^    2 

u  +v 

1  2<J1 

fx(u,v)       =      —^  e  /  t4-1) 

2ttct, 

2^    2 
x    +y 

1  2a2 

f2(x,y)       =      -^2    e  >  <4-2> 

2tt(J2 

2  2 

(u-x)    +(v-y) 

2a2 
D(u-x,v-y)       =      e  .  (4.3) 
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In  (3.35) ,  let 


a 
u 

— 

a 

V 

— 

°i 

a 

X 

= 

a 

y 

= 

°2 

a 

X 

= 

a 

y 

= 

a 

(4.4) 


then  we  have  the  following: 

An  approximate  formula  for  the  optimal  ballistic  disper- 
sion of  the  Simple  Salvo  Model  is  given  by 


2   2        2  2  + 

a0/a   =   [<J>  z   -  z  -  1,  0]  (4.5) 


where 


z   =   (n-l)a^/a2   =   l/(e~*  -  1  +  <j>)  .  (4.6) 


An  upper  bound  to  the  maximum  salvo  kill  probability  is  given 
by 

P  =  1  -  e"(n_1)ni  +  tn-l)(l-e"*)}  .  (4.7) 


2   2  2   2 

The  true  optimal  a^/a      is  a  function  of  n  and  a, /a  ,    but 

2   2 
according  to  (4.5)  and  (4.6)  the  approximate  o2/a      is  deter- 

2   2 
mined  by  z  =  (n-l)a,/a  ,  (j>  being  a  parameter.   To  obtain  an 

2   2 
explicit  formula  for  c~/a   in  z  instead  of  the  parametric 

form,  we  have 
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Lemma  4 . 1 

The  series  solution  to  the  equation  in  <j>, 


e~^    -  1  +  <j>  -  u2/2   =   0  (4.8) 


is  given  by 

<p2      =   u2  +  u  /3  +  u4/12  +  u5/60  +  ...        (4.9) 

Proof 

Let  the  left-hand-side  of  (4.8)  be  expanded  in  <t> . 
Then 


u2      =      4>2(1   -   (j)/ 3   +   <}>2/12    -    <})3/60    +    ...) 


If  the  terms  in  the  right-hand-side  of  (4.9)  are  expanded 
in  (J) ,  we  have 


2       2     3       4        5 
uz   =   <j>*  -  0J/3  +  4)^/12  -  <j>  /60  + 


u3/3      =  $3/3    -    (j>4/6      +    4>5/18    - 


u4/12      =  (J>4/12    -    4>5/18   +    .. 


u5/60      =  c))5/60    -    .  . 


from  which    (4.9)    is   obvious.  D 


56 


Using  the  Lemma  and  the  relationship  (4.6),  we  can  express 

2   2 

a ~ /a   in  (4.5)  as  a  power  series  in  /2/z .   If  we  keep  only 

the  first  three  terms  in  the  expansion,  we  get  an  alterna- 
tive approximate  formula  for  the  optimal  ballistic  disper- 

2 
sion  which  we  will  denote  as  a  . 

a 


aa/ct        =      I   /z/2    "   f  +    3T  /2/z    '  (4.10) 


with 


2      2 
(n-l)a,/a 


In  Table  6  two  values  of  the  approximate  optimal  ballis- 
tic dispersion  are  compared.   One  is  calculated  with  (.4.5) 
and  (4.6),  the  other  with  (4.10).   Agreement  seems  quite 
satisfactory,  differences  appearing  at  most  in  the  fourth 
decimal  digit.   The  formula  (4.10)  needs  more  terms  to  be  a 

Table  6.   Comparison  of  Two  Formulae 


2  ,    2 

2  .    2 

2  ,    2 

z   =    (n-1) 

cr^/a 

a2/a 

a    /a 

800.0 

12.5017 

12.5017 

400.0 

8.5971 

8.5971 

200.0 

5.8367 

5.8367 

100.0 

3.8855 

3.8854 

50.0 

2.5069 

2.5067 

25.0 

1.5335 

1.5331 

12.5 

0.8475 

0.8467 

2   2 
o'/a   is  calculated  by  (4.5),  (4.6),  and 

2   2 
a  /a   is  calculated  by  (4.10)  . 
a7 
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good  approximation  for  smaller  values  of  z,  but  later  we 
will  see  that  the  case  with  z  less  than  25  is  not  important 
from  the  viewpoint  of  application. 

Note 

For  n  =  2,  the  formula  (4.10)  gives 


Qa/a    =   T  °l/a  "  6  +  30  a/al  '  (4'U) 


In  Chapter  II,  the  exact  formula  of  the  optimal  ballistic 

dispersion  was  derived  for  n  =  2.   The  series  expansion  of 

2   2 

a2/a   is 


2   2      1         3    7 
a2/a   =   ^  a1/a  -4  +  54-  a/ai  "  •••         (4.12) 


It  is  interesting  to  note  that  every  term  in  (4.11)  is  smaller 
than  the  corresponding  term  in  (4.12)  :   The  approximate 
value  given  by  (4.11)  for  n  =  2  is  always  smaller  than  the 
true  value. 

For  comparison's  sake,  we  will  also  consider  the  classi- 
cal approximate  formula  by  Merritt  and  King  [Ref .  2]  who 

2 
approximate  the  optimal  ballistic  dispersion  by  a     where 


a2 /a2     -  i  /z~  (4.131 

C  2    c 


z 


=   na^/a   .  (4.14) 
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A.   ACCURACY  INVESTIGATION 

In  this  section  we  will  investigate  the  accuracy  of  both 
our  approximate  formula  and  the  classical  one. 

Table  7A  for  z  =  400,  Table  7B  for  z  =  100,  and  Table 
7C  for  z  =  25  illustrate  the  calculation  carried  out  for  the 

accuracy  study.   For  a  fixed  z,  the  approximate  value  for 

2 
the  optimal  ballistic  dispersion  divided  by  a   is  a  constant 

which  is  given  next  to  the  value  of  z  in  the  table.   On  the 

other  hand,  the  true  optimal  value  increases  with  n:   The 

first  column  is  n  and  in  the  second  column  is  given  the  corres- 

2   2 
ponding  true  optimal  a«/a  ,  the  calculation  of  which  is  based 

on  the  necessary  condition  of  optimality  (2.13),  and  requires 

fairly  long  calculation  time.   The  relative  error  of  our 

2      2    2    2 
approximate  formula  Aa   =  (a  -a0)/a~   x  100%  is  given  in  the 

a      a    z,        <l 

3rd  column.   Its  value  is  found  to  always  be  negative.   The 
salvo  kill  probability  associated  with  o   =  0  is  denoted  as 
P_  and  is  in  the  fourth  column.   The  fifth  column  is  the  upper 
bound  P  of  the  maximum  salvo  kill  probability  given  by  (4.6). 
and  (4.7) .   In  the  sixth  column,  we  have  the  true  maximum 

salvo  kill  probability  designated  simply  as  P  which  is  calcu- 

2 
lated  by  (2.21),  (2.22)  with  the  optimum  a-  given  in  the 

2nd  column.   The  entries  in  the  7th  column  are  the  salvo  kill 

probabilities  associated  with  the  approximately  optimal  dis- 

2 
persion  a  .   The  notation  P   is  used  for  it.   The  last  column 
a  a 

is  the  relative  loss  in  the  kill  probability  which  would  be 

2 
caused  if  we  were  to  use  the  approximate  value  a   instead  of 

a 
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Table  7   Accuracy  of  Approximation  (4.10) 


A 

z  =  400, 

crAl/<**  = 

8.60 

n 

(T*1/*1 

av-^x 

Po 

P 

P 

P<* 

AP% 

5 

9.34 

-7.9 

0.0225 

0.0414 

0.0410 

0.0410 

0.0 

0 

9.46 

-9.1 

0.0634 

0.1484 

0.1456 

0.1454 

0.1 

5 

9.57 

-10.2 

0.1088 

0.2777 

0.2709 

0.2704 

0.2 

0 

9.67 

-11.1 

0.1556 

0.4063 

0.3952 

0.3943 

0.2 

5 

9.76 

-11.9 

0.2024 

0.5229 

0.5081 

0.5067 

0.3 

0 

9.84 

-12.6 

0.2484 

0.6230 

0.6056 

0.6038 

0.3 

5 

9.92 

-13.3 

0.2931 

0.7061 

0.6873 

0.6851 

0.3 

0 

9.98 

-13.9 

0.3362 

0.7733 

0.7542 

0.7518 

0.3 

B 

z  =  100 

,  o^/cx*  = 

3.89 

n 

oV/<*1 

*^x 

Po 

P 

P 

P« 

<*?% 

4 

4.31 

-9.9 

0.0600 

0.0905 

0.0891 

0.0890 

0.1 

8 

4.41 

-11.8 

0.1703 

0.2947 

0.2872 

0.2866 

0.2 

2 

4.48 

-13.4 

0.2828 

0.4952 

0.4809 

0.4794 

0.3 

6 

4.55 

-14.6 

0.3877 

0.6558 

0.6371 

0.6350 

0.3 

0 

4.61 

-15.7 

0.4817 

0.7727 

0.7529 

-  0.7503 

0.3 

4 

4.66 

-16.6 

0.5641 

0.8534 

0.8348 

0.8321 

0.3 

8 

4.70 

-17.4 

0.6352 

0.9071 

0.8910 

0.8885 

0.3 

2 

4.74 

-18.1 

0.6959 

0.9419 

0.9288 

0.9266 

0.2 

n 


z  =  25,  o^oc*  =1.53 


<r2a/<x* 

**?% 

Po 

P 

P 

P«. 

*P/. 

1.77 

-13.6 

0.0573 

0.0660 

0.0652 

0.0651 

0.1 

1.82 

-15.9 

0.2137 

0.2733 

0.2663 

0.2657 

0.2 

1.86 

-17.8 

0.3706 

0.4822 

0.4674 

0.4658 

0.3 

1.90 

-19.3 

0.5076 

0.6497 

0.6299 

0.6274 

0.4 

1.93 

-20.5 

0.6206 

0.7710 

0.7497 

0.7468 

0.4 

1.95 

-21.5 

0.7107 

0.8539 

0.8339 

0.8309 

0.4 

1.98 

-22.4 

0.7812 

0.9085 

0.8912 

0.8884 

0.3 

1.99 

-23.1 

0.8356 

0.9435 

0.9295 

0.9271 

0.3 
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2 
the  true   optimal    a    ,    AP   =    (P  -  p    ) /p  xioo%.      It   is   noted   that 

2.  a 

in  the  case  with  a  large  z  value  the  salvo  kill  probability 
is  much  improved  by  employing  the  optimal  ballistic  disper- 
sion. For  instance,  when  z  =  400,  the  ratio  P/Pn  is  around 
2.5.  But  the  improvement  is  not  so  much  when  z  gets  small. 
In  the  case  with  z  =  25,  the  ratio  is  at  most  1.25. 

As  to  our  approximation,  we  observe  the  following. 

2 

(1)  The  approximate  a   is  always  smaller  than  the  true 

value,  at  least  for  the  values  of  z  investigated. 

(2)  The  discrepancy  between  the  approximate  and  the 

true  values  is  larger  for  smaller  z. 

2 

(3)  For  a  fixed  z,  the  optimal  ballistic  dispersion  a ~ 

increases  as  n  is  increased.   Thus  the  discrepancy 
grows  when  n  increases  for  a  given  z. 

(4)  However,  we  may  say  that  the  discrepancy  is  not  so 
serious  because  the  kill  probability  associated  with 
the  approximate  optimal  ballistic  dispersion  is  not 
so  different  from  the  maximum  value.   The  relative 
loss  given  in  the  last  column  is  at  most  0.4%  which 
occurs  when  z  is  25,  a  less  important  case. 

Figure  5  is  presented  here  to  illustrate  the  properties 

2 
1  and  2.   The  ordinate  is  the  relative  error  Aa  ,  but  it 

a 

2 
needs  several  words.   The  relative  error  Aa   varies  with  n 

a 

for  a  fixed  z  as  seen  in  Table  7.   To  get  a  representative 

2 

Aa   value  for  a  z,  an  n  is  picked  for  the  given  z  such  that 
a 

the  maximum  salvo  kill  probability  is  close  to  0.6,  an 
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2 
arbitrarily  chosen  figure.   Then  the  corresponding  £o   is 

a 

chosen  for  that  z  value,  for  instance  -12.63  for  z  =  400, 
and  similar  procedure  follows  for  other  z's.   The  bottom 
curve  in  Figure  5  shows  this  relationship.   Later  we  will 
discuss  other  curves  in  the  same  figure. 

From  the  figure  we  see  that  the  relative  error  is  always 
negative  and  its  absolute  value  is  decreasing  in  z.   The 
z-dependence  is,  however,  relatively  small,  and  we  learn 
its  bias  to  the  negative  side  is  a  characteristic  of  this 
approximate  formula. 

Similar  investigation  was  carried  out  on  the  classical 

approximate  formula  (4.13).   Table  8A  for  z   =  40  0,  Table  8B 

for  z   =  100,  and  Table  8C  for  z   =25  correspond  to  Table  7, 
c  c 

the  only  difference  being  that  z   is  tabulated,  rather  than 

z.   From  the  tables  we  observe  the  following. 

2   2 

(1)  The  approximate  a  /a   is  larger  than  the  true  value 

except  for  the  cases  of  very  large  z  *s. 

(2)  The  discrepancy  between  the  approximate  and  the  true 
value  is  larger  for  smaller  z  .   This  tendency  is 

similar  to  that  of  our  approximation. 

2 

(3)  For  a  fixed  z  ,  the  optimal  ballistic  dispersion  a9 

increases  as  n  is  increased.   Thus,  the  discrepancy 
gets  smaller  when  n  is  increased  for  a  given  z   less 
than  800. 

(4)  However,  the  discrepancy  is  not  serious  because  the 
kill  probability  associated  with  the  approximate 
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Table  8   Accuracy  of  Approximation  (4.13) 


zc=  400,  o-ci/0<2=  10.00 


9.34 
9.46 
9.57 
9.67 
9.76 
9.84 
9.92 
9.98 


7.1 
5.7 
4.5 
3.4 
2.5 
1.6 
0.9 
0.2 


0.0225 
0.0634 
0.1088 
0.1556 
0.2024 
0.2484 
0.2931 
0.3362 


0.0414 
0.1484 
0.2777 
0.4063 
0.5229 
0.6230 
0.7061 
0.7733 


0.0410 
0.1456 
0.2709 
0.3952 
0.5081 
0.6056 
0.6873 
0.7542 


0.0410 
0.1455 
0.2709 
0.3951 
0.5080 
0.6056 
0.6873 
0.7542 


0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 


B 


zc=    100,  c-cV«2=  5.00 


4.31 
4.41 
4.48 
4.55 
4.61 
4.66 
4.70 
4.74 


15.9 
13.5 


11 
9 

8 
7 
6 
5 


5 
9 
5 
3 
3 
4 


0.0600 
0.1703 
0.2828 
0.3877 
0.4817 
0.5641 
0.6352 
0.6959 


0.0905 
0.2947 
0.4952 
0.6558 
0.7727 
0.8534 
0.9071 
0.9419 


0.0891 
0.2872 
0.4809 
0.6371 
0.7529 
0.8348 
0.8910 
0.9288 


Pc 

0.0889 
0.2866 
0.4800 
0.6363 
0.7523 
0.8344 
0.8907 
0.9286 


*P?e 


3 
0 
9 
2! 
0 


0.1 
0.0 
0.0 


zc=    25,     Ocl/<*1   =    2.50 


<r2a/cc2 

1.77 

1.82 

1.86 

1.90 

1.93 

1.95 

1.98 

1.99 


41.0 
37.1 
34.1 
31.7 
29.7 
28.0 
26.5 
25.3 


0.0573 
0.2137 
0.3706 
0.5076 
0.6206 
0.7107 
0.7812 
0.8356 


0.0660 
0.2733 
0.4822 
0.6497 
0.7710 
0.8539 
0.9085 
0.9435 


0.0652 
0.2663 
0.4674 
0.6299 
0.7497 
0.8339 
0.8912 
0.9295 


0.0648 
0.2638 
0.4631 
0.6249 
0.7451 
0.8300 
0.8881 
0.9272 


AP% 

0.5 
0.9 
0.9 
0.8 
0.6 
0.5 
0.3 
0.2 


64 


optimal  dispersion  is  not  so  different  from  the  true 

maximum.   The  relative  loss  is  at  most  1% ,  and  is 

2   2 
very  small  when  a. /a   is  large. 

The  top  curve  in  Figure  5  is  the  relative  error  Ao   of 
the  classical  approximate  formula,  where  the  abscissa  should 
be  read  as  z   in  this  case.   In  spite  of  its  simple  form,  it 

is  indeed  an  excellent  approximation,  in  particular  for  large 

2 
z  's.   The  z  -dependence  of  the  relative  error  Aa  '  is,  how- 

ever,  more  sharp  than  ours.   The  middle  curve  in  Figure  5  will 

be  discussed  in  the  next  section. 


B.   IMPROVEMENT  OF  THE  APPROXIMATION 

In  the  preceding  section,  we  investigated  the  accuracy 
of  the  approximate  formula  given  by  (4.10) ,  and  it  seemed 
that  (4.10)  gave  a  lower  bound  to  the  true  optimal  ballistic 
dispersion.   The  reason  it  gives  a  possible  lower  bound,  and 
a  method  of  obtaining  improved  approximations  will  be  con- 
sidered now. 

The  Walsh  optimal  f~    and  others  in  (3.16)  are  functions 
of  x, 


x   =  {u*    +    oj2)a2/2  ,  (4.15) 


for  the  Simple  Salvo  MOdel,  and  so  let  us  denote  them  as 
f2(x)  instead  of  f2(oo,,co2),  etc.   With  this  notation,  (3.18; 
is 
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f2(x)   =   p*(x)/D(x)  .  (4.16) 


2   2 
Equation  (3.27)  by  which  the  approximate  a  /a   is  determined 

a 

reads 

2  ,  2 

-xa  /a 

f2(x)   =   e    a     =   pQ(x)/D0(x)  ,  (4.17) 


and 


D(x)   =   DQ(x)  ,  (4.18) 


because  the  damage  function  is  assumed  to  be  Gaussian. 

The  function  p* (x)  is  the  double  Fourier  transform  of 
the  Walsh  optimal  p  function  given  in  Corollary  3.2  with 
a      =  a   =  a-,  ,  and  p.  (x)  is  an  approximate  formula  for  it 
given  by  (3.26) .   Figure  6  illustrates  what  the  curve  of  p* (x) 
is  like.   The  curve  decreases  as  x  is  increased  from  zero, 
changes  its  sign,  and  then  swings  back  to  positive  with  small 
amplitude,  and  so  forth.   On  the  other  hand,  the  approxi- 
mating function  p_ (x)  is  equal  to  p* (x)  at  x  =  0,  has  the 
same  tangent  there,  and  decreases  exponentially.   pn (x)  is 
larger  than  p*(x)  for  a  wide  interval  of  x  as  is  seen  in 
Figure  6 . 


pQ  (x)   >   p*  (x)  ,    0  <  x  <  xQ  . 
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Therefore,  from  (4.16)  and  (4.17),  we  have 


2   2 

-xo  /a  £ 

e   a     >   f2(x)  ,     0  <  x  <  xQ 


2 
which  suggests  that  a   is  smaller  than  the  true  optimal 

value . 

One  of  the  simplest  ways  to  get  an  approximate  formula 

2 
would  be  to  define  a, 

b 

2   2 

-9xa,  /a 

e  pQ(x)/D0(x)  (4.19) 


with  0  <  9  <  1,  or 


a'l      =   a2/0  ,    0  <  6  <  1  ,  (4.20) 

o      a 


and  to  determine  the  parameter  9  empirically  to  achieve 

2 
satisfactory  agreement  between  the  approximate  a,     and  the 

2  2 

true  a~.   Note  that  the  definition  of  a,     corresponds  to 

approximating  the  p* (x)  function  by  p, (x)  with  an  exponent 
1/9  times  the  exponent  of  p« (x) . 

For  example,  let  1/9  be  1.15.   Then, 


at      =   1.15  a2    .  (4.21) 

b  a 


In  Table  9A,  Table  9B,  and  Table  9C,  we  present  the  accuracy 
investigation  of  (4.21).   Accuracy  is  found  much  improved. 
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Table  9    Accuracy  of  Approximation  (4.2/). 


z  =  400,  (TV2/*2  =9.89 


n 

<ra2/c** 

A0-b2% 

Pa 

P 

P 

Pb 

*?% 

5 

9.34 

5.9 

0.0225 

0.0414 

0.0410 

0.0410 

0.0 

10 

9.46 

4.5 

0.0634 

0.1484 

0.1456 

0.1455 

0.0 

15 

9.57 

3.3 

0.1088 

0.2777 

0.2709 

0.2709 

0.0 

20 

9.67 

2.3 

0.1556 

0.4063 

0.3952 

0.3952 

0.0 

25 

9.76 

1.3 

0.2024 

0.5229 

0.5081 

0.5031 

0.0 

30 

9.84 

0.5 

0.2484 

0.6230 

0.6056 

0.6056 

0.0 

35 

9.92 

-0.3 

0.2931 

0.7061 

0.6873 

0.6873 

0.0 

40 

9.98 

-1.0 

0.3362 

0.7733 

0.7542 

0.7542 

0.0 

B         z    =    100 

,    <r*/*z  = 

4.47 

n 

<raVoc2 

*<rba% 

Po 

P 

P 

Pb 

^P/( 

4 

4.31 

3.6 

0.0600 

0.0905 

0.0891 

0.0891 

0.0 

8 

4.41 

1.4 

0.1703 

0.2947 

0.2872 

0.2872 

0.0 

12 

4.48 

-0.4 

0.2828 

0.4952 

0.4809 

0.4809 

0.0 

16 

4.55 

-1.8 

0.3877 

0.6558 

0.6371 

0.6371 

0.0 

20 

4.61 

-3.1 

0.4817 

0.7727 

0.7529 

0.7528 

0.0 

24 

4.66 

-4.1 

0.5641 

0.8534 

0.8348 

0.8347 

0.0 

28 

4.70 

-5.0 

0.6352 

0.9071 

0.8910 

0.8903 

0.0 

32 

4.74 

-5.8 

0.6959 

0.9419 

0.9288 

0.9286 

0.0 

z  =  25,  (rba/o(2s  1-76 


n 

*!*/«* 

4<rb*% 

P0 

P 

P 

Pb 

A??. 

2 

1.77 

-0.6 

0.0573 

0.0660 

0.0652 

0.0652 

0.0 

4 

1.82 

-3.3 

0.2131 

0.2733 

0.2663 

0.2663 

0.0 

6 

1.86 

-5.4 

0.3706 

0.4822 

0.4674 

0.4673 

0.0 

8 

1.90 

-7.1 

0.5076 

0.6497 

0.6299 

0.6296 

0.1 

10 

1.93 

-8.6 

0.6206 

0.7710 

0.7497 

0.7493 

0.1 

12 

1.95 

-9.7 

0.7107 

0.8539 

0.8339 

0.8333 

0.1 

14 

1.98 

-10.8 

0.7812 

0.9085 

0.8912 

0.8906 

0.1 

16 

1.99 

-11.6 

0.8356 

0.9435 

0.9295 

0.9239 

0.1 
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2 
The  maximum  absolute  relative  error  A  a,  reduces  to  12%,  and 

b 

the  relative  loss  in  the  salvo  kill  probability  is,  at  most 

only  0.0  7%.   The  middle  curve  in  Figure  5  is  the  relative 

2 
error  A  a,  . 

b 
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V.   SALVO  MODEL  WITH  NON  GAUSSIAN  DAMAGE  FUNCTION 

In  this  chapter  we  deal  with  a  salvo  kill  model  with 
circular  normal  errors  and  circularly  symmetric  general 
damage  function.   Under  this  general  assumption,  calculation 
of  the  salvo  kill  probability  is  very  difficult,  and  time- 
consuming  computation  is  needed  for  it.   Here  we  propose  an 
approximate  method  for  calculating  the  salvo  kill  proba- 
bility.  By  this  method  computation  is  very  easy  and  yet 
fairly  high  accuracy  is  expected. 

Assumptions  in  this  chapter  are  as  follows:   A  salvo  of 
n  weapons  is  fired.   The  impact  point  error  is  composed  of 
two  parts /  the  random  bias  common  to  all  n  weapons  and  the 
round-to-round  error.   We  adopt  a  Cartesian  coordinate  sys- 
tem (x,y)  the  origin  of  which  coincides  with  the  center  of 
impact  of  the  n  weapons.   The  random  position  of  the  target 
is  denoted  as  (U,V)  with  respect  to  this  coordinate  system, 
and  (U,V)  is  assumed  to  be  circular  normal  centered  at  (.0,0} 
with  variance  a,  : 


fxtu,v) 

= 

1 

2      2 

u   +v 

o    2 

2al 

,      2    C 

(5.1) 


The  impact  points  of  the  ith  weapon  (X.,Y.)  are  independent 
and  identically  distributed  circular  normal  random  variables 
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2 
centered   at    (0,0)    with   variance   a*: 


2^    2 

x   +y 

1  2a2 

f2(x,y)       =      — i-y  e  z       .  (5.2) 

2ttq2 


The  conditional  kill  probability  of  the  target  by  a  single 
weapon  given  that  the  target  is  at  (u,v)  and  the  weapon  impact 
point  is  (x,y)  is  called  as  the  damage  function  D.   In  this 
chapter,  D  is  assumed  to  be  a  function  of  the  miss  distance 


/     2        2 
r  =  V (u-x)  +  (v-y)   only. 


D   =   D(r)  ,     r2   =   (u-x) 2  +  (v-y) 2  .       (5.3) 


Further,  the  lethal  area  is  denoted  as  A,  and  it  is  assumed 
that  there  exist  finite  "moments"  of  all  degree  as  defined 
in  (5.5) , 


2tt  /   D(r)  r  dr   =   A  ,  (5.4) 

0 


CO 

E(rj)   =  =~  j      rj  D(r)  r dr   <   -  ,   j  =  1,2,...   (5.5) 
A   0 

The  conditional  kill  probability  by  a  single  weapon 
given  (U,V)  being  (u,v)  is  given  by  11. 1)/  repeated  here  as 


OO       00 


pCu,v)   =   /    /    D(u-x,v-y) f2 (x,y)  dx  dy  ,        C5.6) 


—00     — OO 
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and  the  salvo  kill  probability  (1.2)  is  renamed  as  (5.7). 


oo   co  n 


=   1  -  /    /   U-p(u,v)}  f1(u/v)  dudv  .        (5.7) 


—CO    — oo 


In  the  Simple  Salvo  Model  in  Chapter  II  with  Gaussian 
assumption  on  the  damage  function,  we  have  a  simple  expres- 
sion for  p(u, v) : 


2^  2 

u  +v 


2       2(a2+ah 
p(u,v)   =   -^ j  e         z       .  (2.6) 

a  +a2 


In  the  general  model  here,  p(u,v)  given  by  (5.6)  is  a  com- 


r~2         2 
plicated  function  in  vu  +  v  ,  and  computation  of  the  kill 

probability  is  difficult.   To  overcome  this  difficulty, 

we  propose  an  approximate  formula  for  p(u,v),  and  then 

present  a  recursive  method  for  computing  the  salvo  kill 

probability  with  this  approximate  p(u,v)  function. 

A.   APPROXIMATE  FORMULA  FOR  p(u,v) 

In  this  section,  an  approximate  formula  for  the  condi- 
tional kill  probability  p(u,v)  is  presented.   As  a  corollary, 
we  will  get  an  approximation  to  the  so-called  Circular 
Coverage  Function. 

Define  the  double  Fourier  transform  of  p(u,v)  as 


0°   °°    i(co,u+w,,v) 
pCoK/uO   =   /    /    e  p(u,v)  dudv  ,        (.5.8) 

— oo    — oo 
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with  similar  notations  being  used  for  D(x,y)  and  f(x,y). 
Then,  from  (5.6),  we  have 


p(co1/co2)       =      D(co1/oj2)  f2  (uj1,u)    )     .  (5.9) 


The  Fourier  transform  of  the  circular  normal  density  (5.2) 
is 


2 

02/  2  _,_  2> 

f2(a)1,u)2)   =   e  *       .  (5.10) 


In  the  Fourier  transform  of  the  damage  function,  let  the 
exponential  function  be  expanded  and  integrated  term  by  term, 


oo        oo        i(a),x+a>2y) 
D(oj1,o)2)       =      /         /         e  D(x,y)   dxdy 

— oo  — oo 


OO  00  oo 

=        I      j        I        jT   i3  (w1x +co2y)  :D(x/y)  dxdy  . 

j  =  0    — oo    — oo   J 


By  symmetry,  the  odd  terms  are  zero,  so 


D(w,,u>2)    =  A      J      72HTX/        /       (co1x+aj2y)^KD(x,y)  dxdy 
k=0  *  — °°      — °° 


=  A   V   |;^}   1  /  {/   (w  cos  9+w0  sin  9) 2kd9 }r2kD(rlr  dr 
k=0  (2k)  l   A  0   0   x       2 
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Therefore, 


D<V«2)/A  -  1  +   I 


T7KTT   2k(2k-2)  ..-2 "E(r   }  (aJi+w2) 


(5.11) 


since 


2tt 


r  2k 

t—  I   (oo,  cos  9  +  u)~  sin  9)   d9   = 
2tt  -*  .    1        2 


(2k-l) (2k-3) • .-1     > 
2k(2k-2) • • -2     '    =  x 


k  =  0 


The  first  few  terms  in  (5.11)  are 


D(u)1,o)2)/A   =   1  -  y  E(r2)  (co2  +  w2)  +  ^jE  (r4  )  (co2  +  w2)  2 


Now  let  us  define 


D  (a),/a>2)   =   Ail  -  j'^i  +  w2^ 


(5.12) 


and  determine  two  constants  b  and  c  such  that  (5.11)  and 
(5.12)  coincide  up  to  the  fourth  power  of  co  .  '  s .  Then  we 
get 


b   =   |  V{E(r2)  }2  -E(r4)/2  , 


(5.13). 


and 
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c   =   -b/2  +  E(r2)/4  . 


If  (5.10)  and  (5.13)  are  substituted  in  (5.9) ,  the  result  is 


S2/  2^  2\ 
b   2    2    — p-taJn+oJ,) 

p  (w1/W  )   =   All  -  j(«£  +u))e  ,  (5.14) 


where 


S2   =   a2  -   b  +  E(r2)/2  (5.15) 


with  b  given  by  (5.13)  . 

Next,  the  inverse  transform  of  p  (u>,  ,u_)  is  calculated: 

a   x  JL 


oo     oo     -i  (U00,  +V0J-)  ~ 

P    (u,v)       =      j  I         /         e  pa^C0l/aJ2)  d(j0idao2    * 

(2tt)    -oo   -°° 


(5.16) 


The  inverse  transform  of  the  first  term  in  (5.14) 


2,2   2,  ,„ 
-S  (uf-HO/2 

Ae     x  A 


is  easily  calculated  as 


2^    2 

u  +v 

A2e 

2ttS 

2S2 

The  second  term  in  (5.14)  is 
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S2  ,  2x  2.  S2   2^  2, 

Ab.  2^  2,  "— ((01+W2)       .  d IA    -r(a3l+aJ2}; 

-  -7-(tt,+&)~)e  =  -b  ~— {Ae  }  , 


and  its  inverse  transform  is  given  by 


2^  2 

u  +v 


K  d ;  A       2S2  , 

-b  j- { =-  e       } 

d(S  )  2ttS 


Therefore  we   have 


2^   2 

u   +v 

p    (U,v)       =      —^- {1    -   -4d  -3L+5L)  }e      2S         .  (5.17) 

a  2ttS^  S2  2S^ 


It  is  noted  that  (5.17)  is  exact  when  the  damage  function 

is  Gaussian  with  parameter  a;  (5.17)  coincides  with  the  exact 

2 
formula  (2.6)  because  in  this  case  A  =  2ira  ,  b  =  0,  and 

c2    2  .  2 

S  =  a-  +  a  . 

The  cookie  cutter  damage  function  differs  substantially 
from  the  Gaussian,  so  consideration  of  the  former  provides 
a  good  test  of  the  robustness  of  (5.17)  .   The  cookie  cutter 
damage  function  is 


.  I1 
( 


if   0  <  r  ^  a 

(5.18) 
otherwise 
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The  probability 


00       00 


p(u,v)   =   /    /    D(u-x,v-y)f  (x,y)  dx  dy 


—  OO     —00 


2       2 

(x-u)  +(y-v) 

j      /    e       z  dx  dy 


27TQ2    2X  2<  2 
x  +y  =a 


is  called  the  Circular  Coverage  Function.   In  the  literature, 

r~2     2 

the  notation  P  (a/a2  ^/a.)  is  used  for  it,  where  r  =  vu  +  v  . 
A  standard  method  of  calculating  the  Circular  Coverage  Func- 
tion is  the  algorithm  given  by  Brennan  and  Reed  [Ref .  10] 
based  on  an  infinite  series  expansion,  but  if  the  accuracy 
required  is  not  too  high,  it  is  convenient  to  use  closed  form 
approximations.   Standard  handbooks  [Refs.  11,12]  recommend 
the  next  formula  for  a/a^  <  1. 

2 
r 

2 
PA(a/a2,r/a2)   =  -—■  e  2E   ,  (5.191 


with 


I2      =  a22    +  a2/4  .  (5.20) 


According  to  (5.17),  the  Circular  Coverage  Function 

(x-r)  +y 

P(a/a0,r/a0)  =  j±-  /   /   e  dx  dy       (5.21) 

2     2     ^tto2    2   2.2 

x  +y  la 
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can  be  approximated  by 


2 

r 

2       .       2    "  "  2 
P_(a/a9/r/a,)   =   -2^-{  1  -  -4- ( 1  -  -£-,)  } e   ^b   ,  (5.22) 

a     2     2       2S2      S2     2S2 


where 


b   =   a2/4/3  ,   S2  =  a2  +a2(l  -  l//3)/4  ,  (5.23) 


since  A  =  TTa  ,  E(r2)  =  a2/2,  and  E(r4)  =  a4/3  for  the 
cookie  cutter  damage  function. 

In  Table  10,  the  Circular  Coverage  Functio  P(R,d)  is 

tabulated  for  R  =  0.1  (0.1)  2.0  and  d=  0.0  (0.5)  2.5.   In 

Table  11,  the  error  AA(R,d)  =  P,(R,d>  -  P(R,d)  is  shown. 

For  small  R,  the  error  is  very  small;  for  instance  the  error 

is  less  than  10   when  R  <  1.0.   But  when  R  is  increased 

beyond  1,  the  error  grows  rapidly. 

In  Table  12,  the  error  A  (R,d)  =  P  (R,dl  -  P(.R,dl  of 

a         a 

our  approximate  formula  is  given.   The  accuracy  is  found  much 

better  than  the  above-mentioned  formula.   The  absolute  error 

-4         <  -2 

for  d  =  0  is  less  than  5  x 10   when  R  =  1.0,  less  than  10 

even  when  R  =  1.6  and  usually  shows  much  higher  accuracy. 


B.   RECURRENCE  FORMULA  FOR  THE  SALVO  KILL  PROBABILITY 

If  the  conditional  kill  probability  p(u,v)  by  a  single 
weapon  is  given  in  the  form  (5.17),  we  can  easily  compute  the 
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Table  10    Circular  Coverage  Function  P(R,d) 


0.00       0.50       1.00       1.50       2.00       2.50 


0.1 

0.0050 

0.0044 

0.0030 

0.0016 

0.0007 

0.0002 

0.2 

0.0198 

0.0175 

0.0121 

0.0065 

0.0027 

0.0009 

0.3 

0.0440 

0.0389 

0.0270 

0.0146 

0.0062 

0.0021 

0.4 

0.0769 

0.0682 

0.0476 

0.0261 

0.0112 

0.0038 

0.5 

0.1175 

0.1045 

0.0735 

0.0408 

0.0179 

0.0062 

0.6 

0.1647 

0.1470 

0.1043 

0.0589 

0.0264 

0.0094 

0.7 

0.2173 

0.1946 

0.1397 

0.0803 

0.0369 

0.0135 

0.8 

0.2739 

0.2463 

0.1790 

0.1050 

0.0495 

0.0187 

0.9 

0.3330 

0.3009 

0.2217 

0.1328 

0.0645 

0.0252 

1.0 

0.3935 

0.3573 

0.2671 

0.1638 

0.0819 

0.0332 

1.1 

0.4539 

0.4144 

0.3146 

0.1976 

0.1019 

0.0428 

1.2 

0.5132 

0.4712 

0.3635 

0.2341 

0.1247 

0.0543 

1.3 

0.5704 

0.5267 

0.4132 

0.2730 

0.1501 

0.0679 

1.4 

0.6247 

0.5802 

0.4628 

0.3138 

0.1784 

0.0838 

1.5 

0.6753 

0.6309 

0.5120 

0.3563 

0.2092 

0.1021 

1.6 

0.7220 

0.6785 

0.5599 

0.3999 

0.2426 

0.1228 

1.7 

0.7643 

0.7224 

0.6062 

0.4442 

0.2784 

0.1462 

1.8 

0.8021 

0.7625 

0.6504 

0.4887 

0.3161 

0.1723 

1.9 

0.8355 

0.7987 

0.6921 

0.5329 

0.3556 

0.2010 

2.0 

0.8647 

0.8309 

0.7310 

0.5763 

0.3965 

0.2321 
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Table    11        Error 


AA  (R,d) 


0.00 


0.50 


1.00 


1.50 


2.00 


2.50 


0.1 

-0.0000 

-0.0000 

-0. 0000 

-0.0000 

-0.0000 

-0.0000 

0.2 

0.0000- 

0.0000 

0.0000 

-0.0000 

-0. 0000 

-0.0000 

0.3 

0.0000 

0.0000 

0.0000 

-0.0000 

-0.0000 

-0.0000 

0.4 

0.0000 

0.0000 

0.0000 

-0.0000 

-0.0000 

-0.0000 

0.5 

0.0001 

0.0001 

0.0000 

-0.0000 

-0.0000 

-0.0000 

0.6 

0.0004 

0.0003 

0.0000 

-0.0001 

-0.0001 

-0.0000 

0.7 

0.0010 

0.0007 

0.0001 

-0.0002 

-0.0002 

-0.0000 

0.8 

0.0020 

0.0014 

0.0003 

-0.0004 

-0.0003 

-0.0001 

0.9 

0.0038 

0.0027 

0.0005 

-0.0007 

-0.0006 

-0.0002 

1.0 

0.0065 

0.0046 

0.0010 

-0.0012 

-0.0011 

-0.0003 

1.1 

0.0106 

0.0076 

0.0018 

-0.0018 

-0.0019 

-0.0007 

1.2 

0.0162 

0.0118 

0.0030 

-0.0026 

-0.0030 

-0.0011 

1.3 

0.0236 

0.0174 

0.0048 

-0.0036 

-0.0045 

-0.0019 

1.4 

0.0330 

0.0246 

0.0074 

-0.0047 

-0.0065 

-0.0030 

1.5 

0.0447 

0.0337 

0.0109 

-0.0058 

-0.0090 

-0.0046 

1.6 

0.0585 

0.0447 

0.0155 

-0.0068 

-0.0121 

-0.0067 

1.7 

0.0746 

0.0578 

0.0213 

-0.0076 

-0.0157 

-0.0095 

1.8 

0.0929 

0.0728 

0.0286 

-0.0080 

-0.0197 

-0.0131 

1.9 

0.1132 

0.0898 

0.0374 

-0.0077 

-0.0240 

-0.0174 

2.0 

0.1353 

0.1086 

0.0478 

-0.0065 

-0.0286 

-0.0225 
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Table    12 


Error 


Aq(R,d) 


0.00 


0.50 


1.00 


1.50 


2.00 


2.50 


0.1 

-0.0000 

-0.0000 

-0.0000 

-0 . 0000 

-0.0000 

-0.0000 

0.2 

-0.0000 

-0.0000 

-0. 0000 

-0.0000 

-0.0000 

-0. 0000 

0.3 

-0.0000 

-0.0000 

-0.0000 

-0. 0000 

-0.0000 

-0.0000 

0.4 

-0.0000 

-0.0000 

-0.0000 

0.0000 

-0.0000 

-0.0000 

0.5 

-0.0000 

-0.0000 

0.0000 

0.0000 

0.0000 

-0.0000 

0.6 

-0.0000 

-0.0000 

0.0000 

0.0000 

0.0000 

-0. 0000 

0.7 

-0.0000 

-0.0000 

0.0000 

0.0000 

0.0000 

-0.0000 

0.8 

-0.0001 

-0.0000 

0.0000 

0.0000 

0.0000 

-0.0000 

0.9 

-0.0001 

-0.0001 

0.0000 

0.0000 

0.0000 

-0.0000 

1.0 

-0.0003 

-0.0002 

0.0000 

0.0001 

0.0000 

-0.0000 

1.1 

-0.0006 

-0.0003 

0.0000 

0.0001 

0.0000 

-0.0000 

1.2 

-0.0011 

-0.0007 

0.0000 

0.0003 

0.0001 

-0.0000 

1.3 

-0.0019 

-0.0012 

0.0000 

0.0005 

0.0002 

-0.0001 

1.4 

-0.0031 

-0.0019 

0.0000 

0.0008 

0.0003 

-0.0001 

1.5 

-0.0049 

-0.0031 

-0.0000 

0.0012 

0.0006 

-0.0001 

1.6 

-0.0075 

-0.0048 

-0.0001 

0.0019 

0.0009 

-0.0002 

1.7 

-0.0110 

-0.0072 

-0.0003 

0.0027 

0.0015 

-0.0002 

1.8 

-0.0157 

-0.0104 

-0.0007 

0.0039 

0.0023 

-0.0002 

1.9 

-0.0218 

-0.0146 

-0.0012 

0.0053 

0.0034 

-0.0001 

2.0 

-0.0294 

-0.0199 

-0.0021 

0.0071 

0.0049 

0.0001 
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salvo  kill  probability  in  a  recursive  way.   We  first  observe 
that  the  salvo  kill  probability 


OO       00 


=   1  -  /    /   (1  -  p(u,v)  }nf1(u,v)  dudv 


—  OO      — OO 


where  f,  (u,v)  is  the  circular  normal  density  (5.1)  and 

2 


2       2 
p(u,v)   =   (y  +  3  -^-y)  e  2S   ,    r2=  u2  +  v2  ,         (5.24) 

2SZ 


can  be  expressed  as 


P   =   1  -  K(n,0)  ,  (5.25) 


where 


K(n,0)       =      /    {1    -    (y  +3t/p)e"t/p}ne"t   dt    ,  (5.26) 

0 


2   2 

and  where  p  =  S  /a, .   This  is  simply  a  matter  of  substitution 

and  the  introduction  of  polar  coordinates.   We  then  have 


Proposition  5.1 
The  function 


K(n,j)   =   /  {1  -  (y  +  8t/p)e-t/p}ne-(^/p+1)tdt        (.5.271 
0 


83 


which  has  two  integer  arguments  and  three  parameters  \i ,    3/ 
and  p  satisfies  the  following  recurrence  relation: 

n  a  1,2,  ...  ;    j  =  0,1,. . .         (5.2  8) 


with 


K(0,j)   =   jT^  /   J  =  0,1,  ...  (5.29) 


Proof 

From  (5.27) , 


K(n,j)  =/  {1-  iv  +  &t/p)e-t/p}n-1e-{i/p  +  1)t   dt 
0 


-  /  {1-  (y+et/p)e-t/p}n-1(.y  +  3t/p)e-t/pe-(^p  +  1)tdt 
0 


The  first  term  in  the  right-hand- side  is  K(n-l,j).   The  second 
terms  is,  when  integrated  by  parts, 
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/  {1  -  (y  + 3t/p)e-t/p}n-1(y  +  St/p) e^e" (^/p+1) fcdt 
0 


-  £{1.  (y  +  3t/p)e-t/p}ne-(^p+1)t| 

0 


-  |(J/P+D  /  (1-  (y  +  3t/p)e-t/p}ne-(J/p+1)t  dt 
0 


-  3/  {l-(y+3t/p)e-t/p}ne-t/pe-(^/p+1)t  dt 
0 


=   £d-y)n  -  ^±£  K(n,j)  -  3K(n-l,j+l)  . 


Therefore 


K(n,j)   =   K(n-l,j)  +  £(l-u)n  -  2±£  K(n,j)  -  3K(n-l,j+l) 


from  which  (5.28)  is  immediate.   From  the  definition  of 
K(0,j),  it  is  clear  that  (5.29)  holds. 

Note  1 

We  use  the  notations  y  and  p  here  in  the  same  way  as 
in  Chapter  II  without  any  subscripts  because  no  confusion 
is  expected.   Their  definitions  are  different,  but  they  are 
just  the  counterparts  of  y  and  p  defined  in  Chapter  II. 

Note  2 

It  is  worthwhile  to  note  that  the  recurrence  relation 
(5.28)  has  a  similar  form  to  (2.17)  .   The  p  function  (5.24) 
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involves  an  extra  parameter  3  which  was  missing  in  the  Simple 
Salvo  Model  in  Chapter  II,  and  therefore,  it  is  necessary 
to  introduce  another  index  j . 
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VI.   COOKIE  CUTTER  SALVO  MODEL 

Consider  that  n  weapons  are  fired  at  a  target  in  a 
salvo.   The  impact  point  error  is  composed  of  two  parts, 
the  random  bais  and  the  dispersion.   We  adopt  a  Cartesian 
coordinate  system  (x,y)  the  origin  of  which  coincides  with 
the  center  of  impact  of  the  n  weapons.   The  random  position 
of  the  target  is  denoted  as  (U,V)  with  respect  to  this 

coordinate  system,  and  (U,V)  is  assumed  to  be  circular  normal 

2 
centered  at  (0,0)  with  variance  a,. 

2^  2 
u  +v 

1        2al 
fx(u,v)   =   -±-j  e      -1   .  (6.1) 

2ttq, 


The  impact  points  of  the  ith  weapon  (X.,Y.)  are  indepen- 
dent and  identically  distributed  circular  normal  random 

2 
variables  centered  at  (0,0)  with  variance  a-. 

2^  2 

_  x  +y 

1       2G2 
f2(x,y)   =   J  e        *  (6.2) 

2ttq2 


The  kill  probability  by  a  single  weapon  conditional  on 
(U,V)  =  (u,v)  and  (X,Y)  =  (x,y)  is  a  function  of  u-x  and 
v-y,  D(u-x,v-y),  and  in  this  chapter  it  is  assumed  that 
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(    1      if    (u-x)      + 
V   0      otherwise    , 


2  .  , ,2  <  _2 


(v-y)"  <  a' 
D(u-x,v-y)   =  {  (6.3) 


with  lethal  area 


2 

A   =   tt  a 


In  the  following,  we  call  this  model  the  Cookie  Cutter 
Salvo  Model. 

The  conditional  probability  that  the  target  is  destroyed 
by  a  single  weapon  given  the  random  position  of  the  target 
being  (u,v)  is 

00       00 

p(u,v)   =   /    /    D(u-x,v-y) f2 (x,y)  dx  dy  , 

— oo    — oo 

and  with  (6.2)  and  (6.3) ,  we  have 


2       2 
(u-x)  +(v-y) 

2a. 


p(u,v)  =  — ^-j        /  /    e  dxdy  =  P(a/a2A2+v2/a2), 

(6.4) 


2lTa2  x2+y2ia2 


where  P(R,d)  is  the  Circular  Coverage  Function. 
The  salvo  kill  probability 


oo     oo 


—00     —00 


Pc   =   1  -  /    /   {1  -p(u,v)  }nf1(u/v)  dudv  , 
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with  (6.1)  and  (6.4)  is  readily  rewritten  as 


Pc  =  1  -  /  {1  -P(a/a2,/2ta1/a  )  }ne"t   dt  .       (6.5) 


We  carried  out  a  number  of  computations  of  (6.5)  with 
various  sets  of  parameters  a,  ,  a„,  a,  and  n  by  the  Simpson 
method.   A  robust  algorithm  for  computing  the  Circular 
Coverage  Function  P(R,d)  for  wide  ranges  of  R  and  d  was  needed 
For  this  purpose,  we  used  the  Brennan-Reed  formula  in  a  modi- 
fied form. 

Brennan-Reed  Formula: 


P(R,d)   =    T   g  k   , 
nio   n  n 


kn/kn-l   =   d  /(2n)  '   n  =  1'2'  *•* 


2 

-dV2 
kQ   =   e       , 


2 

gn  =   gn_x  -  (R2/2)n  e"R  /2/n!  ,   n  =  1,2, 


i     -R2/2 
g0  =  1  -  e 
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Modified  Formula: 


P(R,d)   =   e"d  /2  I       (l-e"R/2C)h   ,  (6.7) 

n^O  n   n 


hn/hn-l   =   d2/(2n^  '   n  «  1>2>     ••■ 


hQ   =   1  , 


Cn   =   Cn-1  +  an'   n  =  lf2f    '  "  (6*8) 


an/an_1   -   R2/(2n)  ,   n  =  1,2,  ... 


C0   "   a0   "  1    ' 


The  modified  formula  is  especially  of  use  in  preventing 
underflow  possible  in  cases  with  large  R  and  d. 

In  their  review  paper,  Eckler  and  Burr  presented  a  short 
table  of  the  salvo  kill  probability  of  the  Cookie  Cutter 
Salvo  Model.   But  the  accuracy  is,  as  they  admitted,  sup- 
posedly low  because  it  is  obtained  by  Monte  Carlo  simula- 
tion.  In  Table  13,  we  give  the  correct  values.   The  figures 
in  parentheses  are  from  the  table  by  Eckler  and  Burr. 

In  the  following  sections,  comparison  of  the  Cookie  Cutter 
Salvo  Model  and  the  Simple  Salvo  Model  is  dealt  with  and 
then  the  optimal  ballistic  dispersion  is  investigated,  where 
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Table  13   Salvo  Kill  Probability  :  Kisi/( Eckler-Burr ) 

a/or<  0.25  0.50  0.75 

0.6 

1.2 

n   =   2  1.8 

2.4 

3.0 


0.221 

0.228 

0.199 

(0.224) 

(0.235) 

(0.201) 

0.573 

0.587 

0.557 

(0.574) 

(0.582) 

(0.563) 

0.834 

0.836 

0.813 

(0.834) 

(0.832) 

(0.814) 

0.954 

0.952 

0.939 

(0.949) 

(0.951) 

(0.942) 

0.991 

0.990 

0.985 

(0.990) 

(0.988) 

(0.985) 

0.50 

0.75 

1.00 

0.231 

0.213 

0.177 

(0.225) 

(0.214) 

(0.187) 

0.525 

0.551 

0.514 

(0.521) 

(0.551) 

(0.516) 

0.736 

0.772 

0.764 

(0.728) 

(0.773) 

(0.763) 

0.872 

0.894 

0.896 

(0.865) 

(0.892) 

(0.893) 

0.946 

0.957 

0.958 

(0.945) 

(0.956) 

(0.957) 

Q/CT, 
0.4 

0.8 

B      n   =    5  1.2 

1.6 
2.0 


Ox/  0", 

<*1(T<                   0.50  0.75  1.00  1.50 

0.2                 0.221  0.210  0.176  0.115 

(0.209)  (0.201)  (0.176)  (0.113) 

0.4                 0.458  0.525  0.504  0.381 

(0.447)  (0.531)  (0.499)  (0.386) 

0.6                 0.605  0.713  0.740  0.648 

(0.595)  (0.716)  (0.740)  (0.640) 

C   n   =    20           0.8                 0.715  0.818  0.861  0.828 

(0.706)  (0.817)  (0.846)  (0.828) 

1.0                 0.801  0.884  0.922  0.923 

(0.792)  (0.882)  (0.924)  (0.924) 

1.2                 0.865  0.928  0.956  0.966 

(0.856)  (0.926)  (0.958)  (0.970) 

1.4                 0.912  0.956  0.975  0.985 

(0.900)  (0.953)  (0.975)  (0.984) 
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the  approximate  formula  presented  in  Chapter  III  and  the 
approximate  method  of  computing  the  salvo  kill  probability- 
developed  in  Chapter  V  are  found  to  be  useful. 

A.   COOKIE  CUTTER  VS.  SIMPLE  SALVO  MODEL 

The  Cookie  Cutter  Salvo  Model  differs  from  the  Simple 
Salvo  Model  of  Chapter  II  only  in  the  damage  function.   In 
Chapter  II,  we  assumed 


2       2 

(u-x)  +(v-y) 

2  2 
D(u-x,v-y)      e       a 


2 

with  lethal  area  2ira  .   For  comparison's  sake,  the  lethal 

area  of  the  two  salvo  models  are  now  set  equal 


a2   =   2a2  .  (6.9) 


The  two  models  share  the  common  parameters  a, ,  a  ,  a,  and  n, 
only  being  different  with  respect  to  the  shape  of  the  damage 
function.   The  quantities  associated  with  the  Cookie  Cutter 
Salvo  Model  will  be  given  suffix  c  in  the  following  when 
necessary. 

It  is  interesting  to  compare  the  salvo  kill  probabilities 
of  the  two  models.   We  begin  with  three  extreme  cases  in  which 
the  salvo  kill  probability  has  closed  form  expressions. 

(1)   The  case  n  =  1. 

The  formula  (6.5)  with  (6.4)  is  easily  integrated  and 
give 

92 


2    2    2 
-a  /(af +0 

P    =   1  -  e      ±         z   .  (6.10) 

c 


The  corresponding  single  shot  kill  probability  under  the 
assumption  of  the  Gaussian  damage  function  is 


P   =   1 5 ±* 5-  .  (6.11) 

1  +  a/(a^_  +  aZ) 


Proposition  6.1 

The  Cookie  Cutter  Salvo  Model  with  n  =  1  gives  a  higher 
single  shot  kill  probability  than  that  of  the  Simple  Salvo 
Model. 


P    >   P  .  (6.12) 

c 


Proof 

Recall  that  e~x  <  l/(l+x)  for  x  >  0.   Let  a  /(c^+c^)  =  x 
in  (6.10)  and  (6.11),  then  (6.12)  is  obvious.        D 


Generally,  in  the  salvo  with  n  >  1,  this  inequality 
does  not  necessarily  hold  as  will  be  shown  in  the  following. 


(2)   The  case  a,  =  0 

Let  a,  ■*  0  in  (6.5)  ,  then  we  have 

2,   2 
-not  /a~ 

P   =   1  -  e      Z    .  (6.13). 

c 
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The  corresponding  formula  in  the  Simple  Salvo  Model  is  given 
by (2 . 2  3) ,  namely 


P   =   1  -  1/(1  +a2/o^)n    .  (6.14) 


(3)   The  case  a-  ■  0 


Let  a2  ■*■    0  in  (6.4)  .   We  get 


2  ,  2 
-a  /a, 

P„  =  1  -  e       .  (6.15) 


The  corresponding  formula  in  the  Simple  Salvo  Model  is  given 
by  (2.24) : 


p      -i  _      n(n-l)  ♦  ♦  -1 {f.    lfi. 

F      x    (n+A)  (n-l+A)  • •  •  (1+A)  '  l0,iDJ 


with 


2/  2 

A   =   a  /a,  . 


Proposition  6.2 

Between  the  salvo  kill  probabilities  of  the  Cookie 
Cutter  Salvo  Model  and  the  Simple  Salvo  Model  with  the  same 
parameters,  the  following  relationships  hold: 


If  o.  =0,  then 


P    >   P  .  (6.17) 

c 
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If  a2  =  0,  then 


> 


2  ,  2 


Pc   <   P      iff   a*/a    J   1/X   ,  (6.18) 


where  An  is  a  positive  solution  to  the  equation 


0 


n 
X    -      I      ln(l  +  X/j)   =   0  .  (6.19) 

j-1 


Proof 

The  proof  of  (6.17)  is  just  the  same  as  that  of  Proposi- 
tion 6.1.   As  to  the  second  assertion,  let 


P   =   1  -  e 


_  ^-gU) 


in  (6.16).   Then 


n  n 

g(A)   =    I      ln(l  +  A/j)  ,      dg/dA   =    £   l/(j+A) 
j-1  j-1 


i  2     2 

Since  g(0)  =0,  dg/dA |     >  1,  d  g/dA   <  0 ,  we  have 

A=0 


g(A)   J   A      iff     A  J  AQ  . 


where  An  is  given  by  (6.19).   The  relation  (6.18)  is 
immediate.  D 
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As  to  the  general  case,  we  must  compare  the  two  kill 

probabilities  numerically.   To  illustrate  the  general  idea, 

Table  14  is  given.   In  Table  14A,  P   for  n  =  4  is  tabulated, 

and  the  corresponding  Simple  Salvo  Model  kill  probability 

P  is  in  Table  14B.   It  is  observed  that  P  ' s  in  the  column 

c 

2   2 

a. /a      =   0  are  larger  than  P's  as  stated  in  Proposition  6.2. 

■   •  2   2 

P  ' s  in  the  first  row  corresponding  to  a~/a   =  0  is  smaller 

2 
than  P's  except  for  the  case  a, /a   =  0,  since  1/XQ  =  0.216 

for  n  =  4 . 

As  a  general  tendency,  P   is  larger  than  P  when  o0    is 

relatively  larger  than  a,.   P   is  smaller  than  P  when  a? 

is  considerably  smaller  than  a,  ,  as  is  seen  in  the  upper 

right  corner  of  the  table. 


B.   OPTIMAL  BALLISTIC  DISPERSION 

We  carried  out  a  golden  search  calculation  to  get  the 

2   2 

optimal  ratio  a~/a   for  a  number  of  cases  with  different 

2   2 
combinations  of  n  and  a, /a  .   The  results  are  tabulated  in 

2    2 
Table  15.   Table  15A  gives  the  optimal  a-  /a   of  the  Cookie 

Cutter  Salvo  Model,  and  Table  15B  gives  the  associated 

salvo  kill  probability,  namely  the  maximum  salvo  kill  proba- 

2   2 
bility  for  a  given  set  of  n  and  a,  /a  .   It  is  noted  that  this 

calculation  is  very  time-consuming.   The  corresponding 

values  for  the  Simple  Salvo  Model  are  tabulated  in  Table  16. 

Before  starting  comparison  of  the  optimal  ballistic 

dispersion  of  the  Cookie  Cutter  Salvo  Model  with  that  of  the 
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Table  14    Salvo  Kill  Probability,  n  =  4 


Cookie  Cutter  Salvo  Model,  P 


c 


<r,1/  <x  * 
cr^/oc2  0.0  2.0  4.0  6.0  8.0  10.0 


0.0 

1.000 

0.393 

0.221 

0.154 

0.118 

0.095 

0.5 

1.000 

0.605 

0.389 

0.285 

0.224 

0.185 

1.0 

0.982 

0.626 

0.425 

0.320 

0.256 

0.214 

1.5 

0.931 

0.615 

0.435 

0.335 

0.272 

0.229 

2.0 

0.865 

0.592 

0.432 

0.339 

0.279 

0.237 

2.5 

0.798 

0.565 

0.423 

0.337 

0.281 

0.240 

3.0 

0.736 

0.538 

0.412 

0.333 

0.279 

0.241 

3.5 

0.681 

0.511 

0.399 

0.326 

0.276 

0.240 

4.0 

0.632 

0.486 

0.385 

0.319 

0.272 

0.238 

B    Simple  Salvo  Model,  P . 


(T, 

l/QCl 

!*/«*» 

0.0 

2.0 

4.0 

6.0 

8.0 

10.0 

0.0 

1.000 

0.594 

0.382 

0.280 

0.221 

0.183 

0.5 

0.988 

0.621 

0.421 

0.317 

0.253 

0.211 

1.0 

0.938 

0.613 

0.431 

0.331 

0.268 

0.225 

1.5 

0.870 

0.589 

0.427 

0.334 

0.274 

0.232 

2.0 

0.802 

0.560 

0.417 

0.331 

0.274 

0.233 

2.5 

0.740 

0.531 

0.404 

0.325 

0.271 

0.233 

3.0 

0.684 

0.503 

0.390 

0.317 

0.267 

0.230 

3.5 

0.634 

0.476 

0.375 

0.309 

0.262 

0.227 

4.0 

0.590 

0.452 

0.361 

0.300 

0.256 

0.223 
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Table   15    Optimal  Ballistic  Dispersion  and 
Maximum  Probability  (  Cookie  Cutter  Salvo  Model  ) 


Optimal  Ratio    o~a&x  f    <*  * 


in 
^.Vcx2      2  4         8         16        32        64 


2 

0.4 

1.0 

1.7 

2.7 

4.1 

6.1 

4 

0.8 

1.6 

2.6 

3.9 

5.8 

8.6 

6 

1.2 

2.1 

3.2 

4.9 

7.1 

10.5 

8 

1.6 

2.5 

3.8 

5.6 

8.3 

12.0 

10 

1.8 

2.9 

4.4 

6.4 

9.3 

13.4 

12 

2.1 

3.3 

4.9 

7.0 

10.1 

14.7 

B        Maximum  Salvo  Kill  Probability 

p 

C   mex% 

Q-^lu3-  2  4         8        16        32 


2 

0.476 

0.626 

0.782 

0.905 

0.973 

0.996 

4 

0.300 

0.435 

0.601 

0.769 

0.899 

C.971 

6 

0.224 

0.339 

0.492 

0.667 

0.826 

e.935 

8 

0.182 

0.281 

0.420 

0.590 

0.762 

C.897 

10 

0.154 

0.241 

0.368 

0.531 

0.708 

C.859 

12 

0.135 

0.213 

0.329 

0.484 

0.661 

C.823 
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Table   16    Optimal  Ballistic  Dispersion  and 
Maxmum  Probability  (  Simple  Salvo  Model  ) 


Optimal  Ratio  cV  /oc2  . 


n 

^•Vot3      2  4         8         16        32        64 


2 

0.5 

0.6 

1.3 

2.2 

3.6 

5.6 

4 

0.5 

1.1 

2.0 

3.4 

5.3 

8.1 

6 

0.5 

1.4 

2.6 

4.3 

6.6 

10.0 

8 

0.7 

1.8 

3.1 

5.0 

7.7 

11.5 

10 

0.9 

2.1 

3.6 

5.7 

8.7 

12.9 

12 

1.0 

2.3 

4.0 

6.3 

9.5 

14.1 

B     Maximum  Salvo  Kill  Probability 

P 


max 


^x/ocz       2  4         8        16        32        64 


2 

0.450 

0.622 

0.780 

0.904 

0.973 

0.996 

4 

0.293 

0.431 

0.599 

0.768 

0.899 

0.971 

6 

0.217 

0.334 

0.489 

0.665 

0.826 

0.935 

8 

0.173 

0.274 

0.416 

0.588 

0.761 

0.896 

10 

0.144 

0.233 

0.363 

0.528 

0.707 

0.858 

12 

0.124 

0.204 

0.323 

0.481 

0.659- 

0.823 
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Simple  Salvo  Model,  we  recall  the  approximate  formula  (3.22) 
given  in  Chapter  III. 


00  oo 


a\      =       [i  /         /         x2{p*(x,y)    -    D(x,y) }dx  dy  , 0 ] +    .  (6.20) 


— oo  — oo 


In  the  case  of  the  Cookie  Cutter  Salvo  Model, 


00      oo 


\  j        j         x2  D(x,y)  dxdy   =   a2/2  , 


A 

—00    — oo 


whereas  in  the  Simple  Salvo  Model, 


00      oo 


1  9 

£  /    /    x   D(x,y)  dxdy   =   a' 


—00    — oo 


Therefore,  (6.20)  states  that  the  optimal  ballistic  dispersion 

of  the  Cookie  Cutter  Salvo  Model  is  larger  than  that  of  the 

2 
Simple  Salvo  Model  approximately  by  a  /2: 


a2  /a2  -  a2/a2   =   0.5  .  (6.211 

2c        2 


Now  we  subtract  entities  in  Table  16A  from  those  in  Table 

2    2    2   2 
15A  and  get  the  difference  a_  /a  -  o^/o.     which  is  tabulated  in 

Table  17A.   The  values  are  positive  with  only  one  exception, 

and  the  predicted  value  0.5  is  observed  in  the  upper  right 
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Table   17    Comparison  of  Cookie  Cutter  Salvo  Model 

and  Simple  Salvo  Model 


cr^/or 


Difference   of  Optimal    Ratio  Clc/c*1"-    <51*  /  erf2" 


n 

8  16  32  64 


2 

-0.1 

0.4 

0.4 

0.5 

0.5 

0.5 

4 

0.3 

0.5 

0.6 

0.5 

0.5 

0.5 

6 

0.7 

0.7 

0.6 

0.6 

0.5 

0.5 

8 

0.9 

0.7 

0.7 

0.6 

0.6 

0.5 

10 

0.9 

0.8 

0.8 

0.7 

0.6 

0.5 

12 

1.1 

1.0 

0.9 

0.7 

0.6 

0.6 

B  Difference   of   Maximum   P's 

n 

^/ct1  2  4  8  16  32  64 

2  0.026  0.004  0.002  0.001  0.000  0.000 


4 

0.007 

0.004 

0.002 

0.001 

0.000 

0.000 

6 

0.007 

0.005 

0.003 

0.002 

0.000 

C.001 

8 

0.009 

0.007 

0.004 

0.002 

0.001 

0.001 

10 

0.010 

0.008 

0.005 

0.003 

0.001 

C.001 

12 

0.011 

0.009 

0.006 

0.003 

0.002 

C.001 
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corner  of  the  table.   In  the  below  left  corner  of  the  table, 

the  figures  are  not  close  to  the  value  0.5,  but  these  cases 

correspond  to  low  salvo  kill  probability. 

In  Table  17B,  the  difference  of  the  maximum  probabilities 

of  each  model  is  tabulated.   The  figures  are  all  positive, 

and  lead  to  the  speculation  that  P      >  P     in  all  cases 

r  c  max     max 

despite  the  fact  that  P  may  be  smaller  than  P  when  a-  is 

c  2 

not  chosen  optimally.  The  author  has  not  yet  succeeded  in 
obtaining  a  proof  of  this  speculation,  but  it  might  be  the 
counterpart  for  n  >  1  of  Proposition  6.1  for  n  =  1. 

In  Chapter  V,  we  developed  an  approximate  method  of 
computing  the  salvo  kill  probability  which  can  be  applied 
to  the  Cookie  Cutter  Salvo  Model.   The  salvo  kill  proba- 
bility is  given  by 

P   =   1  -  K(n,0)  ,  (5.25) 


where 


K(n,j)       =      /    {1    -    (y  +et/p)e~t/p}ne"(j/p+1)tdt  (5.27) 

0 


satisfies   the  recurrence   relation, 


n         n8 


K(n,j)    ^KC-loi^iH)  »_-J^_K(„-lfj+ii,     (5.28) 


n      =      1,2,...;  j=0,l, 
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with 


K(0,j)  '     7+Q    '         3  "  X'2'  ■•■  (5'29 


:+p 


The  parameters  are  given  by  (5.22),  (5.23),  and  (5.24). 


(a2/S2)  (1  -  a2/2/3  S2)  , 


3   =   a4/2/3  S4  ,  (6.22) 


S2   =  o22    +   a2(l  -  l//3)/2  . 


Using  these  formulae  together  with  the  golden  section 
search  calculation,  we  obtained  Table  18.   Table  18A  gives 
the  approximate  optimal  ballistic  dispersion  of  the  Cookie 
Cutter  Salvo  Model,  and  the  associated  salvo  kill  probability 

is  in  Table  18B. 

2    2 
In  Table  19A,  the  difference  between  the  optimal  a9  /a 

2    2 
given  in  Table  15A  and  the  approximate  a      /a   in  Table  18A 

is  tabulated.   The  difference  decreases  towards  the  upper 
right  corner  of  the  table. 

In  Table  19B,  the  salvo  kill  probability  of  the  Cookie 
Cutter  Salvo  Model  associated  with  the  approximately  optimal 
ballistic  dispersion  mentioned  above  is  given.   It  corres- 
ponds to  the  maximum  P  given  in  Table  15B.   In  spite  of  the 

2    2 
discrepancies  observed  in  Table  19A,  the  approximate  cr   /a 

2         2 
gives  almost  the  same  kill  probability  as  the  true  a9  /a   gives 
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Table   18    Optimal  Ballistic  Dispersion 

and  the  Associated  Salvo  Kill  Probability 
(  Cookie  Cutter  Salvo  Model,  Approximation  ) 


Optimal  Ratio  G\A    /c*1  . 


n 

o-flct3-  2  4        8 


16        32        64 


2 

0.5 

1.0 

1.7 

2.7 

4.1 

6.1 

4 

0.7 

1.5 

2.5 

3.9 

5.8 

8.6 

6 

1.0 

1.9 

3.1 

4.8 

7.1 

10.5 

8 

1.2 

2.2 

3.6 

5.5 

8.2 

12.0 

10 

1.3 

2.5 

4.1 

6.2 

9.2 

13.4 

12 

1.5 

2.8 

4.5 

6.8 

10.0 

14.7 

B     Maximum  Salvo  Kill  Probability 

P« 


n 

v*/*z  2  4         8        16        32        64 


2 

0.475 

0.625 

0.782 

0.905 

0.973 

0.996 

4 

0.297 

0.432 

0.600 

0.768 

0.899 

0.971 

6 

0.218 

0.334 

0.490 

0.665 

0.826 

0.935 

8 

0.173 

0.274 

0.416 

0.588 

0.762 

0.396 

0 

0.144 

0.234 

0.363 

0.528 

0.707 

0.358 

12  0.124  0.204  0.323  0.481  0.659  0.323 
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Table   19  Accuracy  of  the  Approximation 
(  Cookie  Cutter  Salvo  Model  ) 


Error  in  the  Optimal  Ratio  f     cr2C*/o(l  _  <y   2  / 


n 

0*iVo(2      2  4         8         16        32        64 


2 

-0.1 

0.0 

0.0 

0.0 

0.0 

0.0 

4 

0.1 

0.1 

0.1 

0.0 

0.0 

0.0 

6 

0.2 

0.2 

0.1 

0.1 

0.0 

0.0 

8 

0.4 

0.3 

0.2 

0.1 

0.1 

0.0 

10 

0.5 

0.4 

0.3 

0.2 

0.1 

0.0 

12 

0.6 

0.5 

0.4 

0.2 

0.1 

0.0 

B   Salvo  Kill  Probability  Associated  with 
the  Approximately  Optimal  Ballistic  Dispersion. 

n 

<T*/o(a  2  4  8  16  32  64 

2  0.475  0.626  0.782  0.905  0.973  0.996 


4 

0.300 

0.435 

0.601 

0.769 

0.899 

0.971 

6 

0.223 

0.339 

0.492 

0.667 

0.826 

0.935 

8 

0.181 

0.280 

0.420 

0.590 

0.762 

0.397 

0 

0.153 

0.241 

0.368 

0.531 

0.708 

0.359 

2 

0.134 

0.212 

0.328 

0.484 

0.661 

0.323 
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VII.   SALVO  MODELS  WITHOUT  CIRCULAR  SYMMETRY 

The  approximate  formula  for  the  optimal  ballistic  dis- 
persion presented  in  Chapter  III  covers  the  salvo  models 
without  circular  symmetry.   In  this  chapter  the  formulae 
for  computing  the  salvo  kill  probability  of  these  models 
are  dealt  with. 

Suppose  that  a  salvo  of  n  weapons  is  fired  against  a 
target.   The  impact  point  error  is  composed  of  two  parts, 
one  being  common  to  all  the  weapons,  and  another  being 
round-to-round  dispersion.   The  random  target  position  (U,V) 
with  respect  to  the  coordinate  system  (x,y)  is  assumed  to 

be  elliptical  normal  centered  at  (0,0)  with  variances 

2       2 
a   and  a   respectively. 

u      v    *  2 


2 

2 

u 

V 

n      2 

^    2 

1 

2a 

2a 

e 

u 

V 

2tt  a    a„ 

u    V 

fx(u,v)   =   .,„  .  -   e         v  .        (7.1) 


The  impact  points  of  the  ith  weapon  (X.,Y.)  are  indepen- 
dent and  identically  distributed  elliptical  normal  random 

2       2 
variables  centered  at  (0,0)  with  variances  a   and  a  , 

x      y 

respectively. 


- 

2 

X 

„    2 

2 

.  y 

r,        2 

1 

e 

2a 

X 

2a 

y 

2tt 

ax 

ay 

f9(x,y)   =  j-  -±    e    A  *    .        (7.2) 

2  2tt  a  a.. 
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The  kill  probability  D  by  a  single  weapon  conditional  on  the 
target  being  at  (u,v)  and  the  weapon  impact  point  being  (x,y) 
is  a  function  of  u-x  and  v-y  only,  and  is  given  by  the 
elliptical  Gaussian  function  with  parameters  a   and  a  . 
Specifically, 


2  2 

(u-x)  (v-y) 

2a  2a 

D(u-x,v-y)   =   e      x  y  .             (7.3) 


It  might  be  an  appropriate  model,  e.g.,  for  a  ground  target 
vs.  a  weapon  impacting  the  ground  obliquely,  because  in  this 
case  the  scattering  of  the  fragments  and  thus  the  conditional 
kill  probability  D  is  by  no  means  circular  symmetric  as 
assumed  in  the  previous  chapters.   For  brevity,  we  call 
this  model  the  Elliptical  Normal  Salvo  Model. 

The  kill  probability  of  a  single  weapon  conditional  on 
(U,V)  being  (u,v)  is  given  by 


oo      oo 


p(u,v)   =   /    /    D(u-x,v-y) f2 (x,y)  dx  dy  , 


— OO     —00 


and   using    (7.2)    and    (7.3)    we   get 


2  2 

u  v 


2      2  2      2 

ct    a                          2  (a   +a    )  2  (a   +a    ) 

.         .                      uxuy                              xx  y      y                          m    a\ 

p(u,v)    =                       *                   e  J      2                          (7.4) 

Aa2+a2)  (a2+a2) 
xx      y     y 
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The  salvo  kill  probability  is  then  given  by 


2     2 

u    v 


2     2 

00  2a    2a 


P  -  1  -  27T  *  g   /    /   (l-p(u/v)}ne    u    v  dudv  .    (7.5! 


U  v   ~"°°   "*°° 


From  this  equation,  Grubbs  derived  the  following  formula: 

Proposition  7.1  (Grubbs) 

The  kill  probability  of  the  Elliptical  Normal  Salvo 
Model  stated  above  is  given  by 


p  -   A""'1^^    ' 


where 


/  J i    2         2.  ,    2   ,     2, 
=      aa/v(a     +  a    )  (a     +  a    ) 
x  y  x        x        y        y 


/ 


(7.7) 


,    2   ,      2,  ,    2  ,2         2,  ,    2 

p        =      (a     +a    )/a      ,         p        =      (a     +  a    ) /a 
x  x        x        u  Ky  y        y' '    v 


The  Grubbs  formula  is  not  suited  for  calculation  when  n 
is  large,  since  it  is  an  alternating  series.   To  overcome  this 
difficulty,  Breaux  and  Mohler  [Ref .  13]  gave  a  method  for  calcu- 
lating the  kill  probability  based  on  an  expansion  of  (1-z) 
in  Jacobi  polynomials  rather  than  as  a  binomial  series.   The 
series  is  found  to  converge  to  the  true  value  with  less  than 
n  terms,  which  is  very  attractive  for  the  calculation  of 
salvo  kill  probability  with  a  large  n. 
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The  author  has  tried  to  find  alternative  ways  of  calcu- 
lating the  salvo  kill  probability  which  are  effective  even 
for  the  cases  with  large  n's,  but  at  present,  the  obtained 
results  are  not  promising.   The  following  propositions  are 
given  not  as  efficient  algorithms  for  calculating  the  salvo 
kill  probability,  but  as  possible  hints  for  developing  more 
efficient  ways  of  calculation. 

For  the  Elliptical  Normal  Salvo  Model,  we  have  a 
proposition  similar  to  Proposition  2.5.   We  omit  the  proof, 
since  it  closely  parallels  that  of  Proposition  2.5. 

Proposition  7.2 

The  salvo  miss  probability  Q(n)  of  the  Elliptical  Normal 
Salvo  Model  characterized  by  (7.1),  (7.2),  and  (7.3)  with 
the  weapon  number  n , 


Q(n)       =      l-iLJl  /         /       (l-yes      fc    )ne      x  y       ds  dt  (7.8) 

7T  J  J 

— OO  —CO 


satisfies  the  relationship 


Q(n)       =        I       (?)     (l-y)n  k   yk  q(k).    ,  (7.9) 

k=0      K 


where 


2         ^2 


/d    0  oo  oo  2     ,2  — o    s    — o    t 

q(k)       =     Da/        /      a-e"8"*)^     x  y       dsdt    .         (7.10) 


11  -00  -oo 
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In  the  case  of  Simple  Salvo  Model,  q(k)  had  a  simple 
recurrence  relationship,  but  we  have  not  succeeded  in 
finding  such  a  relationship  in  this  case. 

Corresponding  to  the  recurrence  formula  for  the  salvo 
kill  probability  presented  in  Chapter  II,  we  have  the  next 
propositions:   The  series  given  in  Proposition  7.3  is  an 
infinite  series,  but  for  a  small  5 /p  it  will  converge  rapidly 

Proposition  7.3 

The  salvo  miss  probability  of  the  Elliptical  Normal 
Salvo  Model  with  n  weapons  is  given  by 


Q(n)   =  A  -62/p2       ^2k)  *   (6/p)2k  J  (ii,  2k)  ,       (7.11) 

k=0  2  (k!) 


where 


J(n,k)   =   ^  /°°(l-ye"t/p)n  e_t  tk  dt  , 


with 


2    2      2    2 

l.(!!M   +  -2U-*)  , 


p    =   -( — _   +  — T 


a         a 
u         v 


2,2      2^2 
l,ax  +  gx    ay  +  qy. 

5   =   2( —  "     2   }  ' 

au        av 


(7.12) 
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and 


a   a 
x   y 


^(a2  +  a2)  (a2  +  a2) 
x        x        y        y 


Proof 


In    (7.8),    let    s   =   rcosS,    t   =    r  sine,    and    (p    +p    )/2    =    p, 

x      y 

(p    -p    )/2    =    6.       Then 
x      y 


2  2  2 

ps      +    p    t        =      r    (p  +  6  cos  29) 
x  y 


and   therefore  Q(n)    is 


./p    p  °°  2  2      2tt         x    2  _Q 

n/„\  vKx^y  r     ,..  -r    Nn   -pr    rr  -or     cos  29,.-,       , 

Q(n)       =      *-  J     (1  -ye         )    e  {J         e  d9}r   dr 

17  0  0 


=  2j^P^  I      d-ye"r  )ne"pr  IQ(5r2)  r  dr  , 


where  In(z)  is  the  modified  Bessel  function  of  the  first 
kind  with  order  0,  and  has  an  expansion  formula 


oo 

I0(z)   =    I       (^-)2(|)2k  .  (7.131 

u        k=0   *'    * 


Thus 


111 


Q(n)       =      /p    d    /p    /    (l-ye"t/p)ne'tIn(6t/p)dt 
x   Y  0  U 


CO 


(2k)!     fX/    .2k        1  f     ..  -t/p.n   -t.2k 


=    ^P^Pv      2.      T t(<S/p)         ToFTT  i     d-ye      /M)    e      t      dt. 

X   Y   k=0    2K(k!)^  {ZK) '       0 

D 


Proposition  7.4 

The  integral  J(n,k)  which  has  two  integer  arguments  and 
two  positive  parameters  y  and  p, 


CO 

J(n,k)   =  irr  /  (1  -ye"t/p)n  e_t  tk  dt  (7.14) 

*   0 


satisfies  the  following  recurrence  relation. 


J(n,k)   =   -£L-j(n-l,k)  +5£_j(n,k-l)  ,  (7.15) 


n  =  1 ,  2 ,  ...  ;    k  =  1 ,  2 , 


J(n,0)   =   ~-  J(n-1,0)  +  jfij.  (l-y)n  ,   n  =  1,2,  ...   (7.16) 


J(0,k)   »  1  ,   k  =  1,2,  ...        (7.17) 


112 


Proof 


00 


j(n,k)       =       1/     (l-ne-^)11"1    e_t   tk  dt 


1     (    c\  -t/p.n-1    -t/p-t    ,k 

-  y   vT  J     ( 1  -  y  e  )        e     /H        t     dt 

K"      o 


The   first   term  in  the  right-hand-side    is   J(n-l,k) .      If 
the   second  term   is    integrated  by  parts,    then 


J(n,k)       =     J(n-l,k)    -   X  £■    (i-ye"t/p)ne"t   tk| 

k.    n  0 


00 

1     0     r     /i  -t/p. n    -t    ,     .k-1 

+   r-r^J(l-ye      /Me        kt  dt 

k!    n   J  .  M 


00 

1     p    r    ,,  -t/p.n   -t    ,k    ,. 

rp  £  J     (1  -  ye         M)     e        t      dt 


Therefore,  for  k  ^  0,  we  have 


J(n,k)   =  J(n-l,k)  +£j(n,k-l)  -£J(n,k)  . 


From  this  relation  (7.15)  is  immediate.   If  k  =  0, 


J(n,0)   =  J(n-1,0)  +  £  (1  -y)n  -  £  J(n,0)  , 


and  therefore,  (7.16)  results.   The  relation  (7.17)  is  obvious 

□ 
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It  is  noted  that  all  the  terms  in  the  expansion  (7.11) 
are  positive,  and  that  the  integral  J(n,0)  is  the  salvo  miss 
probability  of  the  Simple  Salvo  Model  dealt  with  in  Chapter 
II.   Therefore  we  have 

Corollary  7.5 

A  lower  bound  to  the  salvo  miss  probability  of  the 
Elliptic  Normal  Salvo  Model  characterized  by  (7.1),  (7.2) 
and  (7.3)  with  weapon  number  n  is  given  by 


QE(n)   =  \/l  -  62/p2  Q(n)  ,  (7.18) 

where  Q(n)  is  the  miss  probability  of  the  Simple  Salvo  Model 
by  n  weapons  with 

2    2     2    2 

,  a  +  a    a  +  a 
1 ,    x   x    y  y, 

a        a 
u        v 

2    2     2,2 
n  a  +  a    a  +  a 

*   _  1  /_x x   _y v, 

5   -   j( 2 2 )  ' 

a        a 
u        v 


y  = 


a  a 

x  y 


>/(a2  +a2)  (a2  +  a2) 
x   x    y    y 
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